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RADAR FOR RAINFALL MEASUREMENTS AND STORM TRACKING 2 


Glenn E. Stout! 


ABSTRACT 


It was recognized during World War II that radar is an excellent instru- 
ment for obtaining data on the spatial distribution of precipitation. During 
the past ten years considerable research has been supported by the Depart- 
ment of Defense at various institutions to determine the application of radar 
in Meteorology. Efforts have been concentrated on the quantitative measure- 
ment of precipitation and on the use of radar for the detection and tracking of 
tornadoes, severe winds, and hailstorms. As a result of studies at the 
Massachusetts Institute of Technology, the Illinois State Water Survey, the 
U. S. Weather Bureau, and the Department of Defense, networks of radar are 
10w being installed by various agencies interested in weather. In some cases 
radar sets are being procured by industry to improve short-range forecasts 
of precipitation and weather which are pertinent to their operations. 

Illustrations show that radar occasionally detects the initiation of a new 
Irecipitation area one to six hours before it reaches a regular hourly re- 
20rting weather station. The areal extent of precipitation is fairly accurately 
lisplayed and the relative intensity of precipitation can be determined. 
Shanges in the intensity and nature of the precipitation as well as in exact 
iming of the event are readily determined. Precise quantitative measure- 
nents for all types of rains have not been obtained to date. Despite this 
imitation, the river forecaster, hydrologist, and meteorologist can make 
ise of radar to improve their skill in short-range forecasting of rainfall and 
‘iver -stages especially during flash floods. 

_ By observing severe thunderstorms on radar scopes that are associated 
tornadoes, the meteorologist is learning to track existing storms and 
erfect his short-range forecast of areas for future warning. The loss of life 
luring the 1957 tornado season was greatly reduced through use of radar to 

ge approaching severe storms. 


Note: Discussion open until June 1, 1959. To extend the closing date one month, a 

_ written request must be filed with the Executive Secretary, ASCE. Paper 1901 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 85, No. HY 1, January, 1959. 

5 ‘Presented at the meeting of the Hydraulics Division of the ASCE in eateaee 

_ on February 25, 1958. 

. Head, Meteorology Section, Illinois State Water Survey, Urbana, Ill. 
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INTRODUCTION 


The existing national raingage network does not provide a reliable esti- 
mate of areal rainfall and often hourly or storm data are not available con- 
current with the event to river forecasters, meteorologists, and other users 


of rainfall data. The raingage merely measures the rainfall in the immediate — 


vicinity of the collector. Numerous investigators(1,2) have shown the inade- 
quacy of point measurements. Figure 1 is an example from one of the Illinois 
State Water Survey’s raingage networks which illustrates the variability of 
storm, monthly, and growing season rainfall for a 400-square mile area in 
central Illinois. Records from fifty recording raingages were used in draw- 
ing the detailed isohyetal patterns. 

During World War II it was observed by military personnel that certain 
frequencies of radar are detectors of precipitation. Preliminary studies 
were initiated after the war at the Massachusetts Institute of Technology, 
through support of the U. S. Army Signal Research and Development Labora- 
tory, to determine the capability of radar to detect meteorological parameters. 
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Figure 1. Daily, monthly and seasonal rainfall patterns for 
400-square mile watershed. 
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In 1948 the Illinois State Water Survey became interested in the utility of 
radar for quantitative measurement of summertime precipitation in Illinois. 
Networks of closely spaced raingages were established to define the repre- 
sentativeness of rainfall observations and compare radar-measured rainfall 
with gage-measured point and areal rainfall.(3) 

Numerous studies have been completed in the past eight years to deter- 
mine the utility of radar in measuring rainfall. Preliminary findings of the 
various research organizations are presented in this paper in order to sum- 
marize progress to date. Difficulties have been experienced because re- 
search on quantitative precipitation measurement by radar have had to use 
discarded equipment of improper wavelength, electronically unstable and un- 
satisfactory test equipment. Despite several limitations, many applications 
of radar are apparent. 


Principal of Radar 


Radar is a high frequency radio device. It emits a short intense pulse of 
snergy which is focused into a narrow beam of invisible energy by a rotating 
antenna, much like a search light. This pulse of energy travels at the speed 
of light. If the beam strikes an object, such as an airplane, a building, or a 
Zroup of raindrops in a rain cloud, a small portion of the energy is reflected 
jack as an “echo” to the point of transmission. The return signal is ampli- 
fied and presented on a cathode ray tube, which is referred to as a Plan- 
Position-Indicator. The range and bearing of the object are automatically 
determined. 

Research in this field by previous investigators(4,5) indicated that radar 
2choes from rain clouds are the result of back scattering of radio energy by 
raindrops falling through the atmosphere. The theoretical expression that 
1as been developed indicates that the energy (power) received from raindrops 
may be expressed as follows: 


= Ft iq6 1 
Pr = Kis ENG (1) 


where P, is received power, P; is power transmitted, r is the distance 
‘rom the radar set to the reflecting raindrops, N is the number of raindrops 
er unit volume, d is the diameter of the raindrops, and K absorbs a con- 
stant for the refractive index of water and a number of constants which are 
yarameters of the radar set. 

Radar does not measure rainfall intensity directly, but measures the re- 
lectivity, Z, which is proportional to 2 Nd®. Several expressions for rain- 
all rate in terms of 2 Nd6 have been obtained from raindrop size data. 
Nhen one of these expressions, 


vhere R is rainfall rate in inches per hour and Kg is a constant of propor- 
ionality, is combined with equation (1), the resulting equation written in 
ogarithmic form is, for the APS-15 radar system, 


36 
¥ R=1.4 x 10! eS ; (3) 
=z f e , PE 
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Limitations : 
Before considering any of the research to date, it may be instructive to 

consider some of the assumptions that are made. The theory is based on 

scattering by spherical raindrops whose diameters are much smaller than 

the wavelength of the radar. The absorption of the signal, also known as at- 

tenuation, and the back-scattered energy are both increased by deviations 

from non-spherical drops. Most drops larger than 1.8 mm are non-spherical. 

Another basic assumption is that the radar beam is filled with precipitating 
particles. Because the beam becomes larger as range increases, this as- 
sumption is not always correct for storms at long ranges, since the rain- 
storm may be smaller than the cross section of the beam. Also, since the 
height of the beam under normal refraction conditions is proportional to the 
square of the range, it is evident that at longer ranges the beam is higher in 
elevation and therefore less likely to be filled. At a range of 100 miles the 
height of the center of the beam is approximately 5000 feet. Owing to this 
limitation, present-day radar used for quantitative rainfall research should 
usually be restricted to 100 miles. 

Another consideration of importance in quantitative precipitation measure- 
ments is that of the relationship of Z, radar reflectivity, to R, rainfall rate. 
Based on 1211 minutes of data, collected by the Illinois State Water Survey, 
using a raindrop camera,(6) Table 1 shows how R varies for given values of 
Z under different conditions of drop size. 


TABLE 1 


Z-R Relationships for Various Z 


Z R R 
If Z = 1119R!+52 If Z = 1y9R1-62 
mm® /ram3 mm/hr mmn/hr 
103 933 dee 
10° 19.0 55.5 
10° 86.1 206 


The table shows that the two Z-R relations yield rainfall rates that differ by — 
a factor of 3 for the same measured Z. e 

Another consideration is that radar obtains a far greater Sampling of rain- 
drops in a shorter period than one can collect manually. In two years of 4 
operating the raindrop camera, Jones(6) sampled only 1211 cubic meters of — 


raindrops, while one pulse from the radar at a distance of 30 miles samples ~ 
556 million cubic meters. 
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The fact that the radar is detecting a volume some distance aloft while the 
precipitation and drop size measurements are determined at the surface is a 
problem of indeterminate seriousness. Due to their smaller terminal veloci- 
ty, small drops are carried farther horizontally than larger drops. This 
drift may become serious in a wind field where the vertical wind shear is 
appreciable. Considerable variation of reflectivity and rainfall rate in the 
vertical may occur under vertical wind shear conditions. Evaporation and 
coalescence of the drops may also change the reflectivity and rainfall rate 
between the radar volume and the ground. Thus, an error in surface rain- 
fall estimation can be induced by those physical processes. 

A further consideration is stability of the radar. Normally the receiver 
Sensitivity is subject to fluctuation and should be monitored quite frequently. 

Another rather serious limitation at 3 cm and shorter wavelengths is at- 
tenuation; that is, loss of signal strength due to intervening rain which ab- 
sorbs a part of the transmitted and back-scattered energy. Robertson and 
King(7) determined experimentally that the attenuation is proportional to the 
product of rainfall intensity and distance into the storm. McGill University(8) 
and the Illinois State Water Survey using observed drop size data and Mie 
scattering theory found that attenuation per unit distance was related expon- 
entially to the rainfall rate. Robertson and King’s relation gives higher values 
of attenuation for rainfall rates below 90 mm/hr. Above 90 mm/hr the ex- 
oenential relations give higher attenuation values. 


Projects to Evaluate Radar 


As part of a project at the Illinois State Water Survey, sponsored by the 
U. S. Army Signal Research and Development Laboratory to evaluate the use- 
fulness of radar in estimating quantitative precipitation, areal rainfall 
amounts for one minute were compared with areal mean radar one-minute 
amounts. The rainfall amounts were obtained from a raingage network of 50 
yages in a 100-square mile area located 20 miles northwest of the radar sta- 
‘ion. Radar-indicated amounts were obtained from a calibrated APS-15A, a 
3-cm radar set. Areal means were determined by planimetering the radar 
isoecho contours. 

Figure 2 shows the comparisons between radar power ratio and raingage 
amounts. The data(9) were derived from 95 minutes of summertime rains 
where no echo intervened between the radar site and the network. Cases of 
marked dissimilarity between radar and raingage patterns were also elimi- 
iated from the analysis. A lag up to 2 minutes between radar and raingage 
amounts was allowed to permit the precipitation to fall to the ground from 
yeam height. Considerable scatter in the data is evident. 

From recent analysis of an unusually heavy rainstorm,(10) the following 
‘onclusions regarding quantitative radar observations at 3 cm were pre- 
sented; 


“The CPS-9 radar set is satisfactory for portraying the areal extent of 
ight rainfall, since attenuation effects are limited in light rain. In moderate 
o heavy rainstorms, this set does an excellent job of delineating the forward 
‘dge and the horizontal extent of an approaching storm in the absence of in- 
ervening echoes, as shown in Figure 3 for the October 9, 1954 storm over 
iorthern Illinois. At 1730 CST numerous rain cells are evidenced in areas 
yetween recording raingages. However, at 2030 CST neither the depth 
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Figure 2. Relationship between rainfall intensity 
and radar-received power factor. 


nor the intensity of a storm is accurately portrayed. Consequently, objective 
methods of obtaining accurate quantitative rainfall estimates with 3-cm rada1 
do not appear feasible. However, an experienced radar meteorologist, using 
synoptic weather information and climatological relationships, along with 
knowledge of the limitation of 3-cm wavelength radar, is able to make better 
estimates of the rainfall in an area than can be obtained from analysis of the 
synoptic weather maps alone.” 


In another study the radar change in intensity was compared with the rain: 
gage intensity change. Using 5-minute amounts for cases of insignificant at- 
tenuation, it was found that 88 percent of the time the two parameters had the 
same sign change. 

Using Robertson-King’s equation, the Illinois State Water Survey has cor- 
rected a number of radar measurements of rainfall for attenuation. Cases 
where echo intervened between the radar were eliminated. Also cases of 
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a RADAR ECHOES 1730 C 
150 NAUT. MILES TOTAL RAINFALL 17-1800 C 


/ RADAR ECHOES 2030 C 
/50 50 CMI TOTAL RAINFALL 20-2100 C 


Figure 3. Radar echoes and hourly rainfall amounts for October 9, 
hy 1954 over northern Illinois and Indiana. (solid lines, 
full gain; solid areas, 26-DB reduction. All station 
_ rainfall values are zero unless otherwise indicated.) 
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strikingly dissimilar isohyetalisoecho patterns were eliminated. Results 
based on a comparison of point rainfall indicate no discernible improvement 
in the 3-cm radar’s ability to measure rainfall. Areal mean rainfall compar- 
isons were also made in an effort to reduce observational error, but no sig- 
nificant improvement in measurement error was obtained. The unsuccessful 
efforts with attenuation correction indicate that other factors are important 
besides attenuation. Some of these parameters have been discussed earlier 
in the section on limitations. 

The Weather Radar group at the Massachusetts Institute of Technology, 
under sponsorship of the U. S. Army Signal Research and Development 
Laboratory, has been testing its pulse integrator with a 10-cm radar as a 
method of determining intensity of rainfall. The pulse integrator is an elec- 
tronic signal intensity averaging device. Utilizing the ranging unit of the 
radar, the peak voltage of the return signal at a point or along an azimuth is 
electronically averaged. Gating allows the scope to be divided into numerous 
areas so that rainfall values can be determined. Figure 4 which was prepared 
from a MIT report(11) is an example of the type of information readily availa- 
ble. Note the similarity in the data when the radar areal rate is compared 
with a point raingage rate. 

Correspondence with MIT personnel indicates that satisfactory accuracy 
(approximately 50 percent) is obtainable, provided the relation between signal 
intensity on their 10-cm radar and rainfall rate is obtained from a series of 
measurements rather than from purely theoretical considerations. Also, one 
or two sensitive raingages in the region being covered should be used as a 
check in any storm. Figure 5 shows the results of some recent work at 
MIT.(12) For light rainfall rates of less than 1/2 inch per hour, there is a 
similarity in traces between the radar and a raingage. However, an error of 
25 percent is evident much of the time. 

Texas A & M in cooperation with the U. S. Weather Bureau has been inter- 
ested in use of radar data for flood forecasting in small basins. Their 
studies(13) are concerned with the velocity, size, and orientation of a storm 
as it moves over a basin. These studies, employing data from 10-cm radar, 
point out the well-known fact that a storm moving downstream along a basin 
produces a much sharper hydrograph than one moving up the same basin, al- 
though the actual volume of rain may be the same. They recommend deter- 
mining the “velocity of maximum concentration vector” for the individual 
drainage basins. This is defined as the direction and speed of movement of a 
given storm that will produce the maximum crest, other things being equal. 
Radar scope composite photos are made by photographing the scope at 5- 
minute intervals for an hour and making one print from these 12 negatives. 
Quantitative radar-rainfall data are obtained using hourly rainfall data and 
the radar composite. 

In another study the Illinois State Water Survey(14) has compared the radai 
indicated rainfall over the extended time intervals, with the Weather Bureau 
hourly raingage data. Radar-indicated rainfall amounts were summed for 3 
hours at 135 Weather Bureau precipitation stations in Illinois and Indiana. 
Radar isohyets for the 3-hour period were drawn to fit these data. Figure 6 
shows the results of this analysis. The assignment of values to these radar 
isohyets, based on the 3-hour recorded precipitation amounts, was attempted 
but proved difficult and in many cases impossible. A radar isohyet would 
pass through too wide a range of raingage values, so that any calibration was 
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PULSE INTEGRATOR (mm/hr) = 6 
RAINFALL (mm/hr) = © 


Figure 4. Comparative pulse integrator and rainfall values for 
July 10, 1952 (courtesy of Massachusetts Institute of 
Technology). ; 


meaningless. The radar isohyets are labeled as A, B, andC. At times core 
lisplacement caused difficulty in calibration. 

The University of Miami in cooperation with the U. S. Weather Bureau is 
-oOmparing 2-hour composite photos with observed rainfall values, using an 
3P-IM, a 10-cm radar set. At 10-cm there is the advantage of less attenua- 
ion than at 3-cm. Since the radar scope, also known as the Plan-Position- 
ndicator (PPI), is an intensity modulated scope, the more intense areas over 
he period of integration will produce greater darkening of the film than less 
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Figure 5. Comparison of data from SCR615-B radar and 


Hudson-Jardi raingage (courtesy of Massachu- 


setts Institute of Technology). 
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—— RADAR ISOHYET 


——— RAINFALL ISOHYET 
f°) 25 50 MILES 


SCALE 


Figure 6. Comparison of 3-hour total rainfall and radar data 
for January 5, 1955. 


intense areas. Also, a storm of longer duration will expose the film more 
than one of short duration. After many types had been tested, Kodak- 
Ektalure film was selected owing to its slow speed and wide latitude of 
exposure. 

Many photographic and radar problems are involved in measuring rainfall 
quantitatively by radar. Critical control of the photographic developing 
process, film shrinkage, and halation are photographic problems. At close 
range, the scope becomes saturated. A logarithmic receiver instead of the 
conventional linear receiver is one solution to this problem. The fact that 
the velocity of the sweep in azimuth is proportional to range causes echoes at 
close range to appear brighter than those farther away. This effect is being 
compensated for by a filter, the density of which is inversely proportional to 
range. Furthermore, critical control of radar performance and setting must 
be ed 

re 7 shows an opacity-range-rainfall nomogram prepared from these 
Ba 15) Although the scatter on these diagrams is large, a potentially use- 
ful relation is indicated. This is a rapid and inexpensive way of supplementing 


hydrologic data. 
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various rainstorms (courtesy of the University of 


Figure 7. Comparison between opacity of film and range for 
Florida). 
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Prospects 


Future developments will probably include more use of photographic inte- 
gration. This is a relatively inexpensive and quick way to obtain areal mean 
rainfall over several hours. More use is likely to be made of color presenta- 
tion. Different colors could be used to represent a variety of intensities and, 
thus, enable an operator to judge the intensity distribution at a glance. Stor- 
jage and memory tubes could enable the operator to survey the storm in time 
jand space. Electronic computers may also be used but they are expensive. 

Transmission of an actual scope picture to distant places may ultimately 
replace the coded message of today. Work is already in progress on trans- 
mitting the scope picture by facsimile, employing a low frequency to permit 
use of land lines. Thus, large scale composites can be made at an analysis 
center. 

_ In the radar sets themselves we can expect the development of a greater 
dynamic range, which will eliminate scope saturation from the photographic 
integration process. Circuitry that will compensate for range attenuation 
may be developed. Constant height PPI presentation may add to a better 
understanding of the precipitation mechanism. 

Development of equipment that would make radar information more ac- 
cessible, useful, and reliable is in prospect. Analysis will most likely be 
aimed at determining areal precipitation over a several hour period and 
deveioping models of precipitation patterns. 


Severe Storm Tracking by Radar 


Since 1953 radar has proven itself to be a valuable instrument in detecting 
the areas of severe hailstorms and tornadoes. Incidental to Illinois State 
Water Survey studies of quantitative precipitation using radar, the first com- 
plete radar recording of the development of a tornado causing damage esti- 
mated at $3,000,000 was made.(16) Figure 8 illustrates the PPI for April 9, 
1953, showing a hooked appendage developed on the southwest corner of a 
large thunderstorm. A tornado was centered under the vortex part of the 
appendage. Since then, numerous observations have been recorded by radars 
in Massachusetts, Oklahoma, Texas and Kansas. It appears that the maximum 
detectable range of hooked-shaped echoes is about 50 miles. Actually, the 
parent thunderstorm cloud may be detected at ranges up to 250 miles. 

_ Exact data on the value of radar through better warnings during local 
severe storms is difficult to ascertain, since data on past storms are very 
incomplete. Experience has indicated that radar is very effective in the de- 
tection and tracking of severe storms. The network of radar stations needs 
to be expanded and improved upon. Better communications between the ob- 
server of the radar scope, the forecaster, and the public should be estab- 
lished. Eventually, the loss of life from tornadoes will be negligible. 


CONCLUSIONS 


- Results to date-on the use of radar to measure precipitation should be 
svaluated with respect to the user’s requirements. Research studies indicate 
hat point or areal measurement of precipitation by radar are feasible, but 


f 


- 
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Figure 8. Photos of radar PPI with 10-mile range markers 
showing development of cloud pattern associated 
with the tornado of April 9, 1953. 
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| Subject to discrepancies of considerable magnitude at times. The accuracy 
: seems to be limited by the uncertainty in the empirical relation between sig- 
nal intensity and rainfall rate which are based upon several assumptions. 
: However, one must keep in mind that rainfall amounts from raingages are 
often not representative of areal mean rainfall. 
| The applications of radar in hydrometeorology are numerous. Concurrent 
with the event, the detection, areal extent, relative intensity and changes in 
intensity, duration, storm velocity, and orientation of precipitation with 
respect to the basin are readily recorded for the immediate use of the river 
or the synoptic forecaster. 

Considerably more research and development of equipment should be pur- 
sued to solve some of the problems currently confronting the radar meteorolo- 
gists. 

Radar data will become more valuable as a network of stations is created, 
communication established to handle the enormous amount of data, and tech- 
niques developed to cope with the limitations of the data. 
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INTERIM CONSIDERATION OF THE COLUMBIA RIVER ENTRANCE 


John B. Lockett,!1 M. ASCE 


SYNOPSIS 


To meet.the increased needs of ocean navigation, the Federal project for 
improvement at the mouth of the Columbia River was modified by the River 
and Harbor Act of September 3, 1954, to provide for an entrance channel of 
suitable alignment having a depth of 48 feet at mean lower low water and a 
width of one-half mile. Under the modified project, this channel would be ob- 
tained by dredging and the construction of a spur jetty on the north shore, in 
combination with existing jetties and contraction works, at an additional cost 
currently estimated at $10, 340,000 for new work. The difficulties associated 
with prosecution of this modified project call for a re-evaluation of the natur- 
al forces controlling the regimen of the entrance area. 


INTRODUCTION 


The Corps of Engineers is attempting to secure by vigorous dredging a 
channel at the entrance to the Columbia River with least depth of 48 feet and 
of such width as may be permissible within the limits of available funds and 
dredging plant and the minimum requirements of present navigation. Although 
project depths of only 35 feet are authorized and maintained along the main 
Columbia River channel from the entrance to the mouth of Willamette River 
and thence to Portland, a distance of 117.5 miles, an additional usable depth 
of 13 feet through the entrance is necessary to assure safety to ocean vessels, 
particularly during the heavy sea swells and storms which occur at that lo- 
cality. With this full depth assured, the physical capacity of the entrance to 
carry ocean commerce will be practically unlimited. However, maintenance 
J F 
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of previous project depths and the work of securing the new project di- 
mensions have been made difficult, not only by adverse weather and sea con- 
ditions, but also over the years by the growth of the hazardous Clatsop Spit 
shoals. These shoals, lying along the channel side of the South Jetty, have 
gradually forced the navigation channel toward the North Jetty, as shown in 
Figs. 3 and 4, increasing annual dredging requirements and creating an unde- 
sirable alignment in the sea approach to the Columbia River entrance. De- 
spite these difficulties, a usable channel, although less than full project di- 
mensions, has been maintained. During the period 1946 to 1956, ocean traffic 
through the mouth of Columbia River increased from nearly 7 to over 12 
million tons per year, as shown in the following tabulation: 


YEAR MILLIONS OF TONS 


1946 
1947 
1948 
1949 
1950 
1951 
1952 
1953 
1954 
1955 
1956 


a, one 
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History 


A chart, made by Admiral Vancouver in 1792, clearly indicated the ex- 
istence of a single entrance channel with depth of 4-1/2 fathoms between Cape — 
Disappointment and Point Adams. In 1839, Sir Edward Belcher revealed the 
existence of two channels at the entrance separated by a large shoal area 
which forced the deeper channel northward against Cape Disappointment to 
turn abruptly southward and join a generally shallower channel along Point 
Adams. This shoal area, called the Middle Sands, grew in size until its 18- 
foot depth contour, in 1874, embraced a very large area seaward and south of 
Cape Disappointment, and covered almost all of the area now occupied by the 
present entrance. A portion of the Middle Sands formed an island inside the 
entrance which later separated from the main shoal area and, migrating to 
the north to form Sand Island, restored a single entrance channel to the mouth 
of the Columbia River in 1885. At that time, the main portion of the Middle 
Sands had moved westward to form a southward extension of Cape Disap- 
pour nent and minimum depths of 20 feet prevailed on the ocean bar as shown 
in Fig. 1. 

Pursuant to recommendations of the Board of Engineers for the Permanent 
Improvement of the Mouth of the Columbia River, (1) Congress on August 2, 
1882, adopted the original permanent navigation project for the Columbia 
River entrance. This project provided for a bar channel 30 feet deep to be ob- 
tained by construction of a South Jetty, 4-1/2 miles in length, lying in a gener- 
ally northwesterly direction from its base at Point Adams. Construction of 
the jetty was started in April 1885 and was completed. with four groins along 
its north side, in October 1895, At that time the crest of the structure sloped 
from 12 feet above mean lower low water at its base to plus 10 feet at a point 
1-1/8 miles from the shore, tence to plus 4 feet at its outer end. A total of 
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946,000 tons of rock was placed in the original construction of the South Jetty 
at a total cost of $1,969,000. 

Although shoal areas marking the birth of Clatsop Spit began to form along 
the north side of the structure during its construction, a general improvement 
of navigable depths over the bar was evident in this period. By 1895 mini- 
mum depths of 31 feet were available over the bar on an alignment which, 
under the influence of the jetty, had migrated to the north about three miles. 
These favorable conditions, however, were short-lived because of continued 
migration of the channel to the north; by 1896 near closure of the 30-foot 
depth contour on the bar occurred; by 1898 depths of only 29 feet were avail- 
able; by 1901 entrance depths had deteriorated further; and in 1902, instead of 
one, there were three separate entrance channels, each being about 22 feet 
deep, as shown in Fig. 2. In the meantime, Clatsop Spit shoals had formed a 
hook-shaped island through the South Jetty at its mid-length which later ex- 
tended to form a solid land connection with Point Adams. 

In view of these conditions, and the need for increased depths to provide a 
margin of safety for navigation, Congress, in the River and Harbor Act of 
March 3, 1905, adopted the recommendations of the Board of Engineers on 
Project for Improvement of Mouth of Columbia River.(2) These recommen- 
dations provided for modification of the existing project to secure an entrance 
channel one-half mile wide and 40 feet deep by extension of the South Jetty, 
construction of a North Jetty and groins, and dredging. Extension of the South 
Jetty, completed to a total length of 6.6 miles in August 1913, restored a 
single entrance along its former alignment to the mouth of the river. By 1911 
a depth of 24 feet was again available across the bar. As completed, the South 
Jetty had an average top width of 25 feet and an elevation above mean sea 
level of 10 feet from the shore to the “knuckle” end of the original project, 
from whence the height was increased by stages to an elevation of plus 24 
feet at its outer end. Construction of the North Jetty began in 1913 and by the 
next year, when only about 1 mile in length, entrance depths had increased to 
30 feet. As constructed, the North Jetty had a top width of 25 feet, side slopes 
of 1 on 1-1/2, and crest elevation of 28 to 32 feet above mean lower low water. 
Completion of the jetty to its full length of 2.4 miles in August 1917 was ac- 
companied by further improvement of entrance depths to 40 feet, which depths 

were available through a narrow bar channel. By 1920, the 40-foot depth 
contours had retreated considerably providing an entrance channel about one 
mile in width; and, during the next five years, these contours retreated farther 
providing a minimum entrance depth of 45 feet in 1925. The next two years 
have been recognized as the period when the most satisfactory channel con- 
ditions prevailed at the mouth of the Columbia River. By 1927, minimum 
depths of 47 feet were available over the bar, as shown in Fig. 3. It is inter- 
esting to note that for the first time, depths in excess of 40 feet prevailed 
along the channel side of the North Jetty. 

However, concurrently with improvement of navigation conditions, shoals 
developed behind the North Jetty and formed a solid land mass extending al- 
most to its full length. Clatsop Spit continued to grow until by 1931 its 
westward and northward extension had reduced entrance depths to about 43 
feet. Rehabilitation of the South Jetty was started in 1931 and completed in 
February 1936; then, later in 1936, the outer 500 feet of the structure was 
impregnated with a hot asphaltic mix in an attempt to prevent raveling at this 
point by heavy seas. This expedient was not fully successful and a solid 
concrete terminal was constructed at the outer end in 1941. Further growth 
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of Clatsop Spit forced the navigation channel closer to the North Jetty placing 
deeper water along the channel side of the outer end of that structure. In an 
effort to prevent continued movement of the channel to the north, four perme- 
able dikes were built along the south shore of Sand Island, and Jetty “A”, ex- 
tending southward from Cape Disappointment, was completed in 1939. The 
continued westward movement of Clatsop Spit resulted in the removal of more 
of the shoal area protecting the North Jetty thereby causing some damage 
from undercutting along the channel face of the structure. During 1938 and 
1939 the jetty was rehabilitated and a concrete terminal block was placed at 
its outer end. By 1944, although the 40-foot depth contour opposite Jetty “A” 
had moved southward, this local improvement in conditions was accompanied 
_ by marked growth of Clatsop Spit to the west and north opposite the end of the 
North Jetty. Continued advance of the spit through 1951 coincided with definite 
movement of the 40-foot depth contour along the north side of the entrance 
channel to the east and south and in 1952, for the first time, an inner bar in 
extension of Clatsop Spit had formed between the outer end of the jetties. 

In addition to 410,000 tons of stone used for rebuilding the outer end of the 
original South Jetty, a total of 4,427,000 tons of stone was placed in the ex- 
tension of the South Jetty, completed in 1913, at a total cost of $5,657,000. 
Nearly 3 million tons of stone were used in the construction of the North Jetty 
at a cost of $4,319,000. Reconstruction of the South Jetty involved the 
placement of 2,288,000 tons of stone in the superstructure at a cost of 
$4,930,000. Rehabilitation of the North Jetty required the placement of 
244,000 tons of stone and 234,000 tons were used in the construction of Jetty 
“A” at a total cost of $1,271,000. 

As ample depths prevailed over the ocean bar, dredging under the 40-foot 
entrance project was confined to the face of Clatsop Spit. Amounts dredged 
since 1939 under that project were as follows: 


AMOUNT DREDGED AMOUNT DREDGED 
FISCAL YEAR CUBIC YARDS FISCAL YEAR CUBIC YARDS 
1939 501,300 1949 1,023,800 
1940 1,319,800 1950 927,200 
1941-44 None 1951 1,000,300 
1945 393,100 1952 1,267 ,400 
1946 186 ,900 1953 2,796,800 
1947 483,000 1954 2,141,700 
1948 1,030,200 1955 1,749, 300 


The full project width was obtained only during the last three fiscal years 
listed above. 

Under the provisions of the project adopted in 1954, the Corps of Engineers 
on April 15, 1956 initiated the work of dredging the 48-foot entrance channel; 
and, by the end of the dredging season in October, minimum depths of 48 feet 
over the full channel width of one-half mile through the inner and outer bars 

were obtained. These operations were accompanied by varied amounts of 
apparent scour and shoaling. Heavy shoaling of the channel during the non- 
dredging period between October and April, together with the lack of sufficient 
dredging plant, resulted in the decision to concentrate the 1957 dredging 
‘program to the securing of project depth throughout a channel of 1500-foot 
width along the adopted alignment. These later operations were also ac- 
‘companied by natural scour and shoaling and heavy shoaling again occurred 
‘during the ensuing non-dredging season. Conditions prevailing in June 1957 
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are shown in Fig. 4. The following table shows the results of the 1956 and 
1957 dredging programs and changes occurring within the adopted channel 
since initiation of the work as revealed by periodic condition surveys: 


TABLE 1 


DREDGING PROGRESS, MOUTH OF COLUMBIA RIVER, 48-FOOT PROJECT 


QUANTITY OF QUANTITY OF APPARENT SCOUR OR SHOALING 
MATERTAL IN MATERIAL DURING PERIOD, 1000 C.Y. 
DATE AUTHORIZED DREDGED 
CHANNEL , DURING PERIOD 
L000 C.Y. L000. .C ove SCOUR SHOALING 
Apr 15, 1956 18,911 
1,301 1,428 = 
May 1, 1956 16,182 
1,554 ? 2,746 
May 15, 1956 17,374 
1,371 1,158 “ 
gun Ls 32956 14,845 
1,392 - 1,396 
Jun 15, 1956 14,849 
1,424 “ 274 
Jul 2, 1956 13,699 
1,467 2,242 ¥ 
Jul 19, 1956 9,990 
746 1,403 z 
Aug 1, 1956 7,841 
1,841 32 . 
Aug 26, 1956 5,968 
1,495 é 696 
Sep 15, 1956 5,169 
1,186 a 95 
Oct 5, 1956 4,078 
659 - 3,285 
Nov 7, 1956 6,704 
0 * 1,067 
Jan 14, 1957 7,771 7 
) - 1,728 
Feb 4, 1957 9,499 i 
0 : 2,301 
Mar 4, 1957 11,800 : 
0 257 - 
Apr 6, 1957 11,543 
373 1,459 zs 
Apr 29, 1957 9,711 ; 
766 1,590 - 
Jun 3, 1957 7,355 : 
723 = 
Jul 10, 1957 - 7,022 ec 
755 - 1,219 
Aug 7, 1957 7,486 a 
1,142 " 
Sep 18, 1957 5,115 . po 
150 _ 606 
Oct 10, 1957 5,571 
0 : 
Nov 5, 1957 6,880 1,309 
Dec 11, 1957 8,515 ‘ LE 
0 a 
Jan 21, 1958 9,456 oak. 
Feb 19, 1958 10,811 9 - 1,355 
) 1,575 ™ 


Mar 13, 1958 9,236 
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Statement of the Problem 


During initial construction of the South Jetty, shoal areas, which later 
formed a land connection with Point Adams, gradually developed along the 
outside of the original structure. Almost concurrently, shoaling occurred 
along the inside of the jetty forcing the main channel to the north. This latter 
shoaling intensified with the extension of the South Jetty and with construction 
of the North Jetty so that a large land mass, Clatsop Spit, now lies inside and 
along the South Jetty. This land mass and its surrounding shallow depths 
present a definite hazard to navigation. During the period January 1947 to 
January 1952, the U. S. Coast Guard has reported a total of five vessels 
grounding in the entrance. In addition, owners of tankships carrying petrole- 
um products reported that two of their vessels had cracked in crossing the 
entrance and presumably had touched bottom. Although depths greater than 
48 feet are available in a natural channel inside the outer bar, just northwest 
of the Clatsop Spit shoals, as shown in Fig. 4, this channel is particularly 
hazardous during severe storms and has a greater deflection angle than the 
adopted channel that encroaches on the shallower depths of the spit. It has 
been possible to maintain project depths throughout usable widths during a 
portion of the time along the adopted entrance alignment by vigilant dredging 
when weather and sea conditions permitted, but such work has been costly. 
The problem resolves itself to one of securing the authorized project di- 
mensions along a suitable alignment at costs which will be commensurate with 
the benefits of improvement. 


Factors Relating to the Problem 


The problem of securing dependable depths along a suitable alignment at 
the mouth of Columbia River has presented a challenge to the Corps of Engi- 
neers for many years. The complexities associated with this problem were 
recognized as early as 1882, when the minority report of the Board of Engi- 
neers for the Permanent Improvement of the Mouth of the Columbia River(1) 
indicated that “* * * The peculiar circumstance of heavy sea waves, strong 
and alternating currents, and the prevalence of sand in the vicinity, in unstable 
positions, introduce obstacles of so serious a nature that they may be said to 
place limitations upon the applications of principles at the Columbia River 
** *”, Later, in 1903, the Board of Engineers on Project for Improvement of 
Mouth of Columbia River in its report, which formed the basis for the 40-foot 
entrance project, gave thorough consideration to all phenomena known at that 
time to have a bearing on the achievement of the desired project depths. 
Subsequently, in 1932, an extensive program of current measurements across 
one principal river range, about 5 miles above the outer end of the jetties, in 
combination with other measurements including bottom sampling was under- 
taken by the Corps of Engineers and a report (3) on this program contains 
data on sand sizes, currents, tidal flows, suspended load, and other items. A 
report 4) on the mechanical and petrographic characteristics of material 
forming the river bottom, estuary, and beach at the mouth of the Columbia 
River was prepared by Dr. Edwin T. Hodge, Consulting Geologist, in April 
1934, These reports; studies under taken by Dr. M. P. O’Brien;(5) the Tidal 
Model Laboratory at Berkeley, California; (6) salinity observations by the 
Corps of Engineers in 1936;(7) together with periodic condition surveys by the 
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Corps of Engineers and data concerning tides, winds, temperatures, and re- 
lated phenomena observed by other Federal and State agencies, form the bases 
for present knowledge of forces controlling the regimen of the entrance area. 
From these data there has been extracted the following information pertaining 
to factors which may have a significant bearing on the problem: 


River Stage and Discharge 


The streamflow pattern of the lower Columbia River is influenced by 
spring snowmelt runoff from upland areas and, to a lesser degree, by rainfall 
during the winter months. The low-flow period generally occurs during the 
months of September through November and during this period river stages 
may infrequently be as low as zero, the adopted low water stage at Vancouver. 
At this location, winter rains from December through February usually result 
| in gage readings of 4 to 5 feet. Beginning in March, melting snows gradually 

increase river stages to greater than 20 feet at Vancouver during the latter 
| part of May and through the month of June. River flows then gradually de- 
‘ cline to the low-flow period of September through November. Near the mouth, 
the extreme low river daily flow is estimated at 60,000 cfs, the average flow 
at about 255,000 cfs, the average spring freshet discharge at about 625,000 cfs, 
. and the maximum discharge at about 1,300,000 cfs. Occasionally, winter 
freshets along the lower reach of the river exceed the average spring freshet. 
The average low water slope in the estuary is about 0.05 foot per mile. 


Tides 


Tides at the mouth of the Columbia River have the diurnal inequality typi- 

| cal of the Pacific Coast with a long runout from higher high to lower low 

water. The mean range of tide is 6.5 feet, the range from mean lower low 

| water to mean higher high water is 8.5 feet, and extreme tides range from 

minus 2.6 feet to plus 11.6 feet mean lower low water. Tidal effect extends 

| upriver to Bonneville Dam, 140 miles, during periods of low flow. The tidal 
prism has been estimated at 782 square-mile feet. The current survey of 

1932 revealed the following characteristics of tides at the entrance: 


. 
| 
| 


: 
| TABLE 2 


TIDAL CHARACTERISTICS, MOUTH OF COLUMBIA RIVER 


AVERAGE RIVER AVERAGE DURATION OF AVERAGE TIDAL FLOW 


MONTH DISCHARGE, TIDES IN HOURS IN 1000 CFS 
1932 1000 CFS FLOOD EBB FLOOD EBB 
“April 315 5.43 6.90 1,037 1,345 

May 556 4.95 7.55 971": 1,481 

September 125 byte k 6.44 L257 1,306 


2 EE 
7 ; 


Winds and Storms 
Two weather seasons having wide variance in wind and sea conditions 


generally occur each year at the entrance. The season of severe storms 
usually extends from October to April and storm winds, beginning generally 
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in the south quadrant, gradually veer to the southwest and finally blow them- 
selves out as they move around to the west or northwest. Storms are often 
continuous for several days and are accompanied by extremely heavy seas 
mostly from the southwest. During the summer season, the conditions of 
wind and sea are much less severe than those of the winter season. From 
May through September, winds from the north or northeast begin in the 
morning, increase in force during the day, and diminish at night. Occasional- 
ly, during the winter and spring, there are several days when an easterly 
wind predominates and the weather is cold. 


Waves 


Waves in excess of 40 feet in height have been reported by the Coast Guard 
at the mouth of the Columbia River. Observations of wave conditions between 
June 1934 and May 1935, as reported by Dr. M. P. O’Brien, indicate that the 
predominate wave direction is from the southwest. Weighted percentage of 
observed waves in each direction were: 


DIRECTION PERCENT DIRECTION PERCENT 
North 529 South £7.5 
Northeast 2.0 Southwest 28.1 
East 3.16 West 24.0 
Southeast BYSS) Northwest 15.6 


Sediment 


Although the Columbia River is not considered to be a sediment-bearing 
stream and is relatively clear for the greater portion of the year, it is known 
that some material is carried in suspension during freshets. Measurements 
below the mouth of Willamette River made in 1922 and at Vancouver in 1950 
indicated 120 and 138 parts per million, respectively, in suspension during 
the average freshet. This, however, does not include the bed load of coarser 
material which, available evidence seems to indicate, travels in waves that 
range from 8 to 14 feet high near the mouth of Willamette River but decrease 
in height downstream and have crests from 300 to 800 feet apart. Tests on 
Eureka Bar, midway between the mouth of Willamette River and the sea, re- 
veal that these waves may move from 30 to 50 feet in a single day. 


Littoral Forces 


The significance of the part littoral forces contribute to the problem may 
be evidenced by the tremendous build-up of material behind the North J etty 
which began with initiation of construction in 1913 and kept pace with the jetty 
as it was pushed seaward. About 48 million cubic yards of material have 
been deposited in that area and, except for the part along the beach below high 
tide, is now stabilized by a natural growth of grass and weeds. The fact that 
there has been no similar build-up of materials behind the outer portions of 
the South Jetty appears to point to a predominate littoral drift to the south 
along the coast in the vicinity of the mouth of the Columbia River. 


Currents 


Although no instruments were specifically employed for the purpose of de- 
termining subsurface current directions, the program of current measurement 
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undertaken in 1932 indicated that mean ebb tide velocities were greater than 
those of the flood tides for all stages of the upper river, and during freshet 
stages were much stronger. The highest velocity observed was 12.28 feet per 
second on the ebb tide at one-tenth depth. These observations appeared to 
indicate that at high river stages, ebb bottom velocities were much higher 
than flood bottom velocities for all percentages of time. At intermediate 
river stage, bottom velocities were lower than at high stage with ebb bottom 
velocities perdominating for about half the measurements. At low river 
stage, flood bottom velocities were generally predominant. Considering the 
magnitude of ebb bottom velocities as compared to flood bottom velocities 
and their greater transporting power, it appeared evident at that time that the 
ebb flow was sufficient to maintain a channel between the jetties of at least 40 
feet deep at mean lower low water and that no further contraction of the 
entrance appeared necessary for satisfactory channel maintenance. The 
following is a summary of mean velocities measured in 1932 opposite Clatsop 
Spit at strength of tides: 


Time of Measurement Average Velocity, Ft. per Sec. 
At Low Stage (September) 


Ebbs (Seven tides) 4.85 
Floods (Seven tides) 4.51 


At Intermediate Stage (April) 


Ebbs (Eight tides) 5.06 
Floods (Eight tides) 3.85 


At High Stage (May) 


Ebbs (Seventeen tides) 5.58 
Floods (Seventeen tides) 3.36 


During August and September 1954, the U. S. Navy Hydrograph Office 
undertook a series of observations at the entrance including use of the Roberts 
Radio Current meter to measure current velocities and direction at depths 
near the surface, at mid-depth, and near the bottom. Measurements were 
taken at a total of seven stations, four located along the entrance channel at 
least one mile outside the jetties, one about two miles south of the outer end 
of the South Jetty, and two just inside the jetties. Although these measure- 
ments were made under difficulties and, accordingly, are incomplete, it is 
interesting to note that, while measurements at depths near the surface and 
at mid-depth indicated a general movement of water into the estuary on flood 
tides and out of the estuary on ebb tides, measurements at three and possibly 
four stations appeared to indicate at the lower levels a general movement of 
water into the estuary even on ebb tides. 

} In January 1958, further observations of current direction and velocity, 
ity, and {properatire were made by the U. S. Navy Hydrograph Office in 
and adjacent to the Columbia River entrance. Measurements made at a point 
‘about four miles upstream from the outer ends of the jetties during a river 
discharge of 175,000 cfs indicated strong density currents. Upstream veloci- 
es at the bottom ranging from 0.6 to 4.4 feet per second were observed dur- 
ing times of downstream and slack currents at the surface. 
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Salinity 

Measurements by the Corps of Engineers(7) in 1936 of the specific gravity 
of surface, mid-depth, and bottom samples at 21 stations along the lower 22 
miles of the Columbia River estuary indicated a definite difference in salinity 
between surface and bottom waters. During periods of flood tide on April 9 
and 10, when the river flow was about 150,000 cfs, the maximum observed 
density of surface water was 1.018 as compared to 1.023 at mid-depth and 
1.024 along the bottom. During periods of half to almost full tide on August 5 
and 6, when the river flow was about 165,000 cfs, the maximum observed 
densities were 1.007, 1.020, and 1.023, respectively. These observations ap- 
peared to indicate that on these dates the upstream limits of salinity pene- 
trations were about 15 and 11 miles above the outer ends of the jetties, re- 
spectively. 

A number of salinity measurements were made by the Oregon State College 
and the Oregon Fish Commission(8) on July 25, 1955 along the navigation 
channel in lower estuary area with findings as shown in the following tabu- 
lation: 


TABLE 3 


SALINITY OBSERVATIONS, COLUMBIA RIVER ESTUARY, .1955 


STA 
(MILES ABOVE MOUTH) SALINITY 0/00 (PARTS PER THOUSAND) AND DEPTH 

AND TIME SURFACE _15' 20° 25° 35' 45" BOTTOM 
STA 1 (2.6 MI) 

0930 1.9 = ae 4.8 “<5 B16 2505 

1535 4.8 = 14.5 == 29.0° se 33.3 
STA 2 (4.8 MI) 

1000 1.5 Ded seme ee ee 24.1 

1600 5.5 2088 fas << 35.8 32.5 
STA 3 (6.8 MI) 

1030 0.8 1.2 5.0 24.3 

1620 2.8 14.8 20.3 23.3 
STA 4 (9.5 MI) 

1100 0.2 0.2 Ot * DE: as 

1625 0.0 0.8 (is 14.8 
STA 5 (12 MI) 

1125 0.0 0.0 0.0 0.0 

1730 0.6 0.8 0.8 3 
STA 6 (15 MZ) 

1200 0.0 0.0 0.0 0.0 

1755 0.4 0.2 0.0 0.0 


Lower high tide occurred at 0528, higher low tide at 1138, and higher high 


tide at 1755 hours on July 25, 1955. River flow was 364,000 cfs on that date. 
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On the following day, these agencies conducted salinity observations at one 
| station in mid-channel off Point Adams Coast Guard Station, Mile 7.5, with 
results as follows: 


TABLE 4 


SALINITY OBSERVATIONS NEAR POINT ADAMS COAST GUARD STATION, 1955 


ES a a a ee ee ee 
TIME OF OBSERVATION AND SALINITY 0/00 


DEPTH 0700 0800 0930 1020 1215 1315 1415 1500 1616 1705 1745 
SURFACE oe. 0.9 0.0 0.0 0.4 0.6 0.3 thsi 1.7 2.4 3.4 
5! -- D.cee O20 ROO OOM TOs 0.5 «elton 255 4 60) 5-304 
10' eee. oet Ol meso Osaem OSS 1.8 45O> 44 40 
15" -- tS sO MERU s 5 se On Ose OS4 2 088 5.56, | eTGr. 36.0 
20! 9.8 3.2 0.2 Lod 0.8 0.4 2.9 Zoi Dis 5.8 Sie 
25° -- Lie3 O29 916.0) 216.1 8.0 5.4 8.2 thee = isn 7p wae) 
| 30° Ef ete OMe te TLS, FeO) 1TS0. 8.0) 17-015. 2 
| 35" Ss Pet see Oot See OOO. 0 al sain 10,0 016,00.17.0 
40' -- TOR Omo nt | 20c8 Me2OCO Oso 1 Oe2) 13951) 14,0 - 1664" 19e0 
| BOTTOM LOSS A195 767.9 = == 20.6 <-- List V4.2" 14.8 1856 2193 
| On that day, July 26, 1955, lower low tide occurred at 0041, lower high tide 
at 0626, higher low tide at 1223, and higher high tide at 1842 hours. River 
| flow was 356,000 cfs. 
| Temperature 
| 


Observations of temperatures at different depths were made by the Oregon 
State College and the Oregon Fish Commission off Point Adams Coast Guard 
| Station on July 26, 1955 concurrently with some of the salinity observations 
i shown in Table 4 as indicated in the following tabulation: 
| 


TABLE 5 


j 
| 
| 


TEMPERATURE OBSERVATIONS NEAR POINT ADAMS COAST GUARD STATION, 1955 


a, 
TIME OF OBSERVATION & TEMPERATURE (CENT) 05 


wd 


0700 0800 1215 1415 _ 
0 

SURFACE 19.6 16.9 i755 17.4 Wis 
5! 16.0 15.9 17.5 17.5 16.9 
10! 15.9 16.0 17.5 ‘ws 16.5 
15! Se 16.0 17.6 173 16.3 
20! 14.2 15.5 fad 17.0 16.2 
' 25! == 13.0 15.0 16.4 14.9 
30! 13.0 12,8 14.4 16.1 14.8 

30 A 

35! = 12.8 14.0. 15.8 14.7 
40! ae 12.9 14.0 14.8 VAY 
_ BOTTOM 12.6 12.9 14.0 14.8 14.7 


q 
. 
' 
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Shoal Materials 


Samples of bottom materials taken at 118 points in the entrance area dur- 
ing the 1932 current survey indicated that the coarser bottom deposits were 
found near the center of the section off Clatsop Spit and that seaward from 
this section the deposits became finer, the finest materials being found on the 
outer slope of the bar beyond the mouth. Bottom samples graded from coarse 
to fine from the navigation channel toward the jetties inside and from fine to 
coarse toward the jetties on the outside. 

As the result of a study based on analysis of 199 samples taken in the vi- 
cinity of the entrance in 1933 and 1934, Dr. Edwin T. Hodge, Consulting Ge- 
ologist, was of the opinion that the adjacent coast has in recent geological 
time been submerged and that the sea is awash against a cliff cut at a much 
earlier date. He believed that the materials derived from the erosion of the 
coast lie on the ocean floor as a submerged coastal plain and that the beaches 
are constructional, built from sand that has come from the upper reaches of 
the Columbia River and not from the adjacent cliffs. The clearness or 
freshness of the materials, their angularity, and the presence of several 
particles adhering together were considered as proof of the transportation of 
the grains by floating during a recent period of time. He concluded that the 
beach sands are not derived from the erosion of bedrock adjacent to the coast 
but originate from a distant source. As the source was considered one of 
about one-fourth andesite and two-thirds granite and rhyolite rock, it was be- 
lieved that the Cascade Gorge could supply the former, but only the upper 
reaches of the Columbia River could supply the latter materials. 

Sieve analyses of bin samples of materials dredged at the entrance within 
recent years show the following characteristics: 


CLATSOP SPIT BIN SAMPLES 


DATE OF SAMPLING DRAIN SIZE, MM UNIFORMITY COEFFICIENT, 60 

July 23, 1952 -20 1.40 Me 

May 15, 1953 -20 1.30 

June 8, 1953 ol 1.29 

June 10, 1954 a2 1.30 

July 26, 1955 ook 1.43 

May 31, 1956 18 1.55 

August 22, 1956 Le 1.18 

August 22, 1956 puts! 2.00 

September 16, 1957 bys 1.42 

OUTER BAR BIN SAMPLES 

DATE OF SAMPLING Dj pGRAIN SIZE, MM UNIFORMITY COEFFICIENT, 260 
10 

July 26, 1955 eZ 1.50 

July 26, 1955 14 1.64 

May 31, 1956 “gil 1.83 

August 22, 1956 alee 2.08 

August 22, 1956 25 1.60 

September 16, 1957 oL3 153 


Near the end of the 1956 dreding season an area of hard packed clay, ap- 
proximately 3000 feet long and 600 feet wide, was encountered at the 48-foot_ 
level along the centerline of the channel near the turn in alignment. 
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Channel Alignment in Estuary 


Referring to Fig. 5, the adopted alignment of the navigation channel down- 
stream from Harrington Point makes a long curve to the south past Tongue 
Point and Astoria, thence through a sharp reverse curve past Clatsop Spit 
shoals to enter the Pacific Ocean through the jetties in a general south- 
westerly direction. Although this alignment is dictated in part by harbor fa- 
cilities at Astoria and vicinity, it appears possible that a change in alignment 
so as to direct the major river flow due west from Harrington Point past 

) Point Ellice and thence to the entrance might be a factor which could bene- 
ficially influence dredging requirements of the entrance channel and could 
also provide for greater ease in navigation. Such a modification in alignment 
would, of course, require the maintenance of a channel connecting the re- 
aligned main navigation channel with the Astoria Harbor area and adjacent 
facilities. 


Preliminary Analysis of the Problem 


As preface to further consideration of phenomena occurring at the entrance 

‘to Columbia River, review of the aforementioned difficulties to navigation 
| points to the conclusion that there is a definite problem at this estuary 

entrance worthy of detailed study in the interest of orderly and timely plan- 
Ning. In such problems, in which complexities are apparently compounded by 
unknown interrelationships of vaguely defined forces, preliminary analysis, 
based upon available knowledge and data, frequently serves to lend sound di- 
rection to the course of further study. 

Past condition surveys represent the source of much of the available infor- 
mation concerning conditions at the mouth of the Columbia River. Examination 
of these surveys reveals that greater than project depths have prevailed since 
the beginning of improvement throughout a natural channel, with length rang- 

: ing from 4 to 7 miles, lying just inside the bar between Clatsop Spit and Sand 
‘Island. Attainment of project depths across the bar in 1917 provided a direct 
: connection between this deep channel and the sea, a connection which pre- 
vailed until the formation of the inner bar between the jetties in 1952. Con- 

currently, however, the enlargement of Clatsop Spit pushed this channel 
against the North Jetty and made it necessary to dredge the face of the spit in 
order to maintain authorized project dimensions. As indicated previously, 
this dredging, confined to the Clatsop Spit area, totaled nearly 15,000,000 

cubic yards during the period 1939-55. To obtain the increase in project depth 

to 48 feet, dredging since the spring of 1956 has been necessary on both the 

outer and inner bars. 

Comparisons of computed quantities of material remaining within the 
authorized channel, shown in Table 1, with reported bin measurements of 

quantities dredged indicate that recent dredging operations have sometimes 
been accompanied by a natural scouring action which has greatly assisted the 
work of the dredges, and at other times by a heavy shoaling within the channel 
against which the dredges have made little or no headway. Although there ap- 
pears to be no consistent relationship of scour or shoaling to river discharge, 
there does appear to be a tendency for the entrance to scour during a rising 
river condition and to shoal on a falling or almost stationary condition, Dur- 
ing the 150-day period, between October 6, 1956 and March 4, 1957, when the 
river flow was quite uniform and averaged about 140,000 cfs, the entrance 
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hannel shoaled approximately 8,380,000 cubic yards as shown on Table 1. 
his was partially compensated during the next 90-day period to June 3, 1957, 
hen natural scour amounted to nearly 3,300,000 cubic yards. This shoaling 
rend was also evident during the winter of 1957-58 when the amount of ma- 
erial in the authorized channel was increased by a total of 5,850,000 cubic 
oe in the 154-day period between September 18, 1957 and February 19, 
58. 
The major portion of the materials dredged at the entrance during 1956 
nd 1957 has been dumped by the hopper dredges in an area lying about one 
ile south of the extreme end of the South Jetty. Condition surveys reveal 
at, although some 14 million cubic yards of material have been disposed in 
hat location since the spring of 1956, very little of this material has re- 
ined in that area. 
Another phenomenon noted, particularly by condition surveys of April 1957, 
d later in November of that year, was a formation of waves of material, 
nging from 3 to 10 feet in height, along the bottom of the entrance channel 
n both the inner and outer bars. The degree of wave formation within the 
ventrance channel area appeared more pronounced than in any of the adjacent 
festuary bottom areas shown by these condition surveys. 
As pointed out above, sieve analyses of bin samples of dredged materials 
indicate that these materials have the general fineness characteristic of sand, 
though on one occasion material coming within the classification of clay 
as encountered at the 48-foot depth level. Previous analyses of the compo- 
ition of these sands, revealing the presence of andesite, granite, and rhyo- 
ite, led to an early assumption that their source lies in the upper reaches of 
he Columbia River. Although these upper reaches may have been the source 
f a portion of the entrance sands during a previous geological period, the 
legree of shoaling occurring during the winters of 1956-57 and 1957-58 noted 
. when the upland discharge was relatively low, clearly indicates that 
ie upper Columbia River reaches do not now constitute the immediate source 
r locality from which material is now moved into the channel during, say, a 
Single year. A possible source of these materials may be in the littoral drift 
offshore along the coast; the evidence provided by the massive build-up of ma- 
terials behind the North Jetty indicates that these littoral movements, pre- 
dominately in a southward direction, may be of some significance particularly 
when influenced by wave action. In addition to littoral drift, it is also possible 
that offshore bottom areas may be a source of the shoaling materials. An- 
other possibility, although seemingly less plausible, is that these materials 
may have been deposited in the entrance area by wind action. Complete engi- 
neering consideration of the problem should, of course, definitely ascertain 
the immediate source of the shoaling materials. 
_ In considering the possible forces which may be responsible for the 
movement of materials within an estuary area, available knowledge limits 
these forces to river currents, density currents, littoral drift, wind, and 
waves, or any possible combinations of these forces. River currents have 
long been recognized as the primary agent responsible for movement of shoal- 
ing materials in nontidal streams and in the nontidal portions of tidal streams. 
However, the action of river currents is known to be modified considerably 
rithin an estuary by density currents created by differences in specific gravity 
#f the fresh upland discharge and the salt water of the sea. 
mA paper by H. B. Simmons(9) points out that engineers have been keenly 
are in recent years that the amount of fresh water discharged into an 
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estuary and the degree to which it mixes with the salt water of the sea are 
major factors establishing the hydraulic and shoaling regimens of.estuaries. 
With the source of fresh water at the upstream end of an estuary and the 
source of heavier salt water at the seaward end, each type of water tends to 
assume, as the result of density differences, a rough wedge shape with the 
base of the wedge at the source. The sea water entering an estuary does so 
in this general wedge form, the entering tip of which moves along the bottom. 
The interface or line of demarcation between the salt and fresh water may 
vary from well defined to almost obscure depending on the degree of mixing 
of the salt and fresh waters. The ratio of mean upland discharge over a mean 
tidal cycle interval to the mean tidal prism, as measured by the volume of 
water flowing into an estuary from the sea during an average flood-tide peri- 
od, may be used to roughly define three broad categories of mixing: highly 
stratified, when the ratio is 1.0 or greater; partly mixed, when the ratio is 
in the order of 0.25; and well mixed, when the ratio is appreciably less than 
0.1. 

Distribution of current velocities in a highly stratified estuary, such as 
the Southwest Pass of the Mississippi River, shows that upstream from the 
tip of the salt-water wedge the direction of flow throughout the entire depth is 
toward the sea at all times with a vertical velocity distribution similar to that 
of an upland stream. Farther downstream the flow in the fresh-water strata 
is still toward the sea, but the direction of flow in the underlying salt water 
is upstream at all times. The cause of the continuous upstream flow within 
the salt-water wedge is constant erosion of salt water from the interface by 
the outflowing fresh water. The salt water thus eroded must be replaced to 
maintain stability and, accordingly, a continuous upstream current is generat 
ed beneath the interface. In this type of estuary, rapid shoaling usually oc- 
curs in the region of the salt-water intrusion tip. The heavier particles of 
material, which are rolled or pushed along the river bed by the outgoing frest 
water, come to rest as soon as the tip of intrusion is reached. The lighter 
particles, transported largely in suspension in the fresh water, gradually fall 
through the interface along the length of intrusion and are transported up- 
stream by salt water currents to the vicinity of the tip. Flocculation of dis- 
solved materials in the fresh water may also occur along the interface and 
form part of the shoaling mass. 

Available records indicate that upland discharge into the Columbia River 
estuary has ranged from an extreme low of 60,000 cfs to a maximum of about 
1,300,000 cfs. The ratios of these discharges to measurements of average 
flood-tide flows shown in Table 2 indicate that the mixing of fresh and salt 
water within this estuary may assume characteristics of the highly-stratified 
estuary during the spring freshets and of the partly-mixed estuary during 
normal and low-flow periods. Evidence that there exists pronounced salinity 
differences between the top and bottom layers of water of the estuary is avail 
able from observations made by Oregon State College and Oregon Fish Com- 
mission, shown in Tables 3 and 4, on July 25 and 26, 1955, at which time the 
upland discharge was about 360,000 cfs. This finding was confirmed by obser 
vations of the U. S. Hydrograph Office made in 1954 and 1958 as noted above. 
Earlier observations in 1936 indicated that the upstream penetration of the 
salt-water wedge may range from 11 to 15 miles. 

As the result of a recent preliminary reanalysis of the findings of the 
Columbia River current survey of 1932, Simmons, in a later paper( 10) stated 
that the effects of density are of major importance with respect to shoaling at 
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e Columbia River entrance. While recognizing that more detailed studies 
ater may show that certain of the resulting indications are not completely 
alid, this preliminary reanalysis, made for three different upland discharges, 
revealed that at a point 5 miles above the outer ends of the jetties, the flow 
{ 90 per cent depth was predominantly upstream during an upland discharge 
f 127,000 cfs and essentially in balance (50 per cent upstream and 50 per 
ent downstream) during upland discharges of 320,000 and 584,000 cfs. Con- 
sideration of only those current velocities in excess of 0.6 feet per second, 
competent for movement of sand with median grain size of 0.4 to 0.5 mm (bin 
Sample analyses noted above indicate median grain sizes of 0.25 mm on the 
inner bar and 0.20 mm on the outer bar) resulted in a marked increase in the 
upstream predominance of effective flow at 90 per cent depth for an upland 
discharge of 127,000 cfs, a change from a balanced flow condition to one of 
upstream predominance for a discharge of 320,000 cfs, and no change from a 
balanced flow condition for a discharge of 584,000 cfs. Consideration of only 
those current velocities in excess of 1.5 feet per second, at and above which 
velocity the rate of sand movement should be appreciable, resulted in a 
jfurther increase in flood predominance at 90 per cent depth for upland dis- 
charges of 127,000 and 320,000 cfs, and a change from a balanced flow con- 
idition to one of definite upstream predominance for a discharge of 584,000 cfs. 
The above available evidence and reanalysis of the 1932 current survey, to- 
gether with information not susceptible to complete check at this time, point 
to great probability that density currents, created by salinity differences, 
constitute one of the major, if not the most dominant, forces controlling the 
shoaling regimen of the Columbia River entrance. 
Except for its possible effect in creating the large build-up of materials 
|behind the North Jetty, very little information is available concerning the 
nature or extent of offshore littoral drift movements at the Columbia entrance. 
\Considering the great mass of materials now trapped and consolidated behind 
that structure, it may reasonably be assumed that an appreciable amount of 
material is moved by littoral forces across the entrance area. It may also 
be reasonable to expect that a portion of this material may be carried through 
the entrance by flood tides and waves and deposited so as not to be removed 
‘by ebb tides. Clatsop Spit and Sand Island might have had their origins in 
such action and study to determine the mineral and chemical composition and 
particle-shape characteristics, as well as mechanical analyses, should be 
made of the materials forming these geographic features to ascertain whether 
or not they are formed of littoral materials. 
_ Although the existing project provides for the construction of spur Jetty 
“B” to reduce maintenance costs for dredging and to provide protection 
against damage to the shorelines and main jetties from scour which may be 
caused by a change in currents resulting from the deeper channel, prelimi- 
nary analysis of the shoaling problem raises doubt concerning the effective- 
ness of that structure, particularly when viewed in the light of the probable 
rces controlling the shoaling regimen of the entrance area. Considering the 
s jetty structure as a contraction work whose greatest effectiveness would 
occur during the high-flow periods, it is doubted that Jetty “B” would be 
instrumental in reducing, or affecting to any material degree, heavy shoaling 
such as occurred during the relatively low-flow periods in the winters of 
1956- -57 and 1957-58. These factors, as well as the cost of the structure, 
urrently estimated at $5,944,000, and the possible economies which might be 
ptiected by alteration of its design in the light of a detailed appraisal of the 
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forces to be encountered, point to the advisability of a more considered ap- 
proach to solution of the problem. 

Preliminary analysis indicates that such an approach should include a 
comprehensive study of all phenomena susceptible of contributing to shoaling 
within the Columbia River entrance area, of which the phenomena of density 
currents and littoral drift appear to be the primary and controlling agents. 
Such a study should be directed toward answering the following questions: 


a. What and where are the immediate sources of shoaling materials? 

b. What are the effects of density currents and littoral drift on shoaling ? 

c. If density currents and littoral drift are primary agents for moving ma- 
terials from the sea into the estuary, what counteracting forces create 
the balance of conditions in the entrance area? 


Plan of Action 


In the light of experience gained in the study and improvement of tidal 
entrances elsewhere in the United States, the problem at the Columbia River 
entrance was presented to the Committee on Tidal Hydraulics of the Corps of 
Engineers early in 1956 and as the result of that presentation the Committee 
recommended that a program of prototype measurements be undertaken as 
soon as possible with the view to investigating the problem in a hydraulic 
model at an early date. 

Available data which might be useful in a model study were examined, and 
based on this examination, the Director, Waterways Experiment Station, 
Vicksburg, Mississippi in March 1957, recommended a program of hydraulic 
and salinity measurements required for the desigi and verification of a hy- 
draulic model. This program consists of three cycles of measurements to 
obtain data for conditions of low, normal, and high river discharge, each cycl 
to include observations of current direction and velocity for a continuous full 
tidal period of about 25 hours at 19 stations located on six ranges as shown 
in Fig. 5. At each station measurements will be made on 30-minute intervals 
at the surface and bottom and at the intervening quarter points of depth. 
Simultaneous observations will be made initially at one station on each range 
for a continuous period of about 25 hours to establish the relationship betweer 
ranges followed by simultaneous observations at all stations on each range 
until all observations have been obtained. Salinities will be observed at hourl 
intervals concurrently with the observations of current velocities and di- 
rections. If practicable, an attempt will be made to obtain a few measure- 
ments of current velocities, current directions, and salinities on a range be- 
tween the jetties. Considering the equipment which may be available for this 
purpose, it is anticipated that a complete measurement cycle can be complet- 
ed in six days. 

It is planned to undertake the above program of prototype measurements 
as soon as necessary instrumentation, operating personnel and equipment, an 
funds can be made available. It is recognized that, even with an expedited 
program of prototype measurements followed by the immediate design, con- ; 
struction, and verification of a hydraulic model, it may be a number of years 
before the forces that control the regimen of the Columbia River entrance 
may be accurately defined and suitable corrective measures designed and con 
structed. In the meantime, it is anticipated that sufficient maintenance funds 
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will be made available to provide, by dredging, existing project depths over 
usable widths throughout the entrance channel. 


CONCLUSION 


An engineering approach to solution of the problem at the Columbia River 
entrance should include a qualitative definition of forces controlling the regi- 
men of the estuary area, such as may be obtained by means of a prototype 
measurement program, together with a quantitative analysis of the magnitudes 
and effects of these forces as may be provided by means of a hydraulic model 
study. Such steps in the past have led to the solution of complex tidal 
problems in the Delaware River, Charleston Harbor, Savannah Harbor, and 
other improved river estuaries by revealing the structural measures neces- 
sary to make effective use of the forces controlling the particular estuarine 
regimens involved and similar steps should be no less instrumental in point- 
ing to a solution of the Columbia River entrance problem. Even though a 
model investigation might conclusively show that solution of this problem lies 
not in additional structural measures but in regular channel dredging, this 
| investigation would be of inestimable value in that it would definitely elimi- 
| nate further structural measures from future consideration. In the light of 

such a finding, substantial savings could be effected by limiting future expendi- 
tures to maintenance dredging operations only. 


—— 
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AIR MODEL STUDIES OF HYDRAULIC DOWNPULL ON LARGE GATES2 
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| ABSTRACT 


Air model studies were used to determine the hydraulic downpull forces 
on a large fixed-wheel gate and on a slide gate. The applicability of air tests 
Hor this work is discussed, and the method of applying the data is detailed. 
Data are presented that will be helpful in evaluating downpull forces on other 
|= gate structures. 
| 


INTRODUCTION 


The use of air models for studying hydraulic problems has repeatedly 
proved to be a reliable and accurate expedient. The advantages of air tests 
jin place of hydraulic tests have been thoroughly discussed by Hunter Rouse(1) 
and others, and many examples of the use of these air models are given by 

J. W. Ball and D. W. Appel.(2) Another unusual use of air models is pre- 
sented in this paper, wherein studies of hydraulic downpull forces ona eat 
|fixed- -wheel gate and on a high-pressure slide gate are discussed. 
The downpull forces being considered here are those forces produced 

when flow occurs beneath the gate leaf, thereby reducing the pressures on 
ithe bottom of the leaf relative to the pressures on the top. This pressure 
difference, acting on the cross-sectional area of the leaf, produces a down- 
pull force which must be considered in addition to frictional forces and the 

leaf weight in designing the stem and hoist. The method used for obtaining 

e downpull was to determine the unit pressures acting on the top and bot- 
‘om surfaces of the gate leaves, and applying these pressures to appropriate 
areas to compute the downpull forces. No measurements of stem loads were 


made f 


Note: Discussion open until June 1, 1959. To extend the closing date one month, a 
7 written request must be filed with the Executive Secretary, ASCE. Paper 1903 is 
_ part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 85, No. HY 1, January, 1959. 
. Presented at the June 1958 ASCE Convention in Portland, Oregon. 
. Hydr. Research Engr., Hydr. Structures and Equipment Section, Hydr. 
Lab., Bureau of Reclamation, Denver, Colo. 
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When the downpull studies were first considered, it was taken for granted 
that model studies would be made with water. But it was later recognized 
that air model studies would be equally satisfactory and that they would be 
less expensive and time consuming. The method of testing with air would be 
identical to that for water because the air velocities would be maintained be- 
low about 250 feet per second, and incompressible flow equations could be 
used.(2) The degree of error would not exceed about 5 percent at a velocity 
of 250 fps and would decrease rapidly with decreases in velocity. The 
velocities actually used were 64 fps or less, and the error was about 1 per- 
cent. Wetting would not be a problem in an air model and, therefore, un- 
treated wood and metal surfaces could be used. The construction could be 
relatively lightweight because the density of air is low and the weight and 
pressure forces are small. The supply and exhaust systems posed no prob- 
lems because the atmosphere served as both reservoir and receiver. 

The argument in favor of hydraulic models was that by their use free 
water surfaces could be obtained. With air models discharging into the at- 
mosphere no free surfaces are possible. With either of the two gates under 
consideration free surfaces will exist during free discharge operation. 

This limitation was not as great a handicap for the air models as was first 
thought. In the case of the fixed-wheel gate, the maximum downpull occurs 
when there is no free surface and where back pressure or submerged condi- 
tions exist. Thus, only tests with submerged flow were required. These are 
accurately reproduced in an air model. 

In the case of the slide gate, the downpull forces were desired at both free 
discharge and submerged conditions. After study of the leaf shape, it was 
found that the pressure acting on the leaf could be analyzed for either free 
discharge or submerged conditions by proper treatment of the data obtained 
from an air model. The treatment is discussed in detail later in this paper. 


Gate Structures 


The first gate studied was a fixed-wheel gate designed for the power con- 
duit at Glendo Dam, Wyoming (Figure 1). This gate has a leaf 3.3 feet thick, 
16.5 feet wide, and 21.0 feet high. It is located near the entrance of the 21- 
foot-diameter outlet and power tunnel and normally remains fully open. 

When a regular closure is required, it is made with no flow taking place and 
with the pressures upstream and downstream from the leaf balanced. Emer- 
gency closures with unbalanced pressures may also be made. Sustained 
operation with unbalanced pressures will occur during the tunnel filling perio 
when the gate will be used as the flow regulator at openings of 3 to 6 inches. 

The control point on this fixed-wheel gate is at the upstream bottom edge 
of the leaf. When the gate is controlling the flow, the pressures on the leaf 
bottom will be less than the pressures on the leaf top and less than in the tun 
nel just downstream. This occurs because as the flow passes beneath the lea 
the flow velocity— and hence the kinetic energy—becomes high. Simultaneous 
ly, and because the total energy in the water remains essentially constant, th 
piezometric pressures become low. 

The pressures on the leaf bottom were of greatest interest during the tun- 
nel filling period when the gate openings were very small. At these small 
openings, the ratio of opening height to gate thickness was small, so that in 
effect, the opening under the leaf approximated a “short tube,” or more 
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specifically, a “short slot.” As the flow passed beneath the leaf, it experi- 
enced it’s minimum cross section at the vena contracta, and then tended to 
expand and occupy most of the section between the conduit floor and the web 
of the bottom beam at the downstream face of the leaf. This action could pre- 
vent aeration between the leaf and the jet during free discharge operation, 

and would restrict the relief afforded by circulation of water during sub- 
merged operation. Thus, lower pressures than usual would occur on the en- 
tire leaf bottom at very small gate openings, and the downpull forces would 
be extremely large. When the gate opening increases, the jet contraction 

and subsequent expansion changes and more space is left between the jet and 
the downstream web. The fluid is readily able to move upstream through the 
space to relieve the low pressures on the gate bottom. When the gate reaches 
the full open position, the pressures beneath the leaf are about the same as 
those in the tunnel and on top of the leaf. At this point, the downpull is quite 
small. 

The second gate studied was the slide gate design now being extensively 
used by the Bureau of Reclamation. This gate is a rectangular, downstream 
seal, high capacity slide gate designed for use as a regulator under high 
heads (Figure 2). The principal features of the gate are the small slots, the 
outwardly offset downstream slot corners followed by gradually converging 
side walls, and a leaf with a flat upstream face and a 45° sloping bottom. This 
gate was developed to meet the requirements of the outlet works at Palisades 
Dam where large flows of water at heads up to 240 feet are controlled. The 
gate has since been included in the designs of several outlet structures with 
heads as high as 373 feet. 

The magnitude of the hydraulic downpull forces on the fixed-wheel gate 
and on the slide gate were desired so that adequate, but not excessive, hoists 
and handling equipment could be provided. Model studies were made on the 
gates to evaluate these downpull forces. 


The Models 


Air was used as the flowing fluid in the models. A centrifugal blower drew 
air from the atmosphere and forced it through a 10-inch pipeline 6.4 diame- 
ters long, and then through the test section and back to the atmosphere 
(Figure 3). The rate of flow was measured by a flat plate orifice at the en- 
trance of the 12-inch-diameter blower inlet line. All of the parts for the 
models were of simple and inexpensive construction, and were fabricated in 
a short time. 

The test section for the fixed-wheel gate consisted of a plywood conduit 
16.67 inches high (Figure 4). The conduit was 8.25 inches wide upstream 
from the gate leaf, and 8.50 inches wide downstream from the leaf. The 2.82- 
inch-thick leaf was made of 12-gage sheet metal and contained 10 piezéme- 
ters on the bottom along or near the conduit center line. Piezometers in the 
tunnel roof upstream and downstream from the leaf were used to obtain the 
pressure drop across the gate. The slots were simplified and represented 
the width but not the depth of the Glendo slots. The position of the gate leaf 
could be observed through transparent plastic windows that formed the outer 
slot walls. 

The test section for the slide gate was 16.67 inches high and 8.25 inches 
wide (Figure 5). The 2.48-inch-thick leaf was made of heavy gage sheet meta 
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Air was drawn from the atmosphere into inlet 
orifice at center by centrifugal blower at left. 
Flow passed through a 10-inch-diameter pipe 
section, then the test section, and back into 
the atmosphere. 


FIGURE 3 - AIR MODEL AND TEST FACILITIES 


and was supported in slots 0.64 inch wide. Three rows of eight piezometers 
were placed on the leaf bottom so the pressure distribution could be deter- 
mined at distances of 0.18 and 1.10 inches from the sides of the leaf, and on 
the center line. An additional piezometer was placed on the bottom of one of 
the projections of the leaf that extends into the gate slots. The upstream 
pressure head was measured by a piezometer in the right wall 19 inches up- 
stream from the leaf. A piezometer on the roof center line 1 inch upstream 
from the leaf gave the approximate bonnet pressure. No bonnet was included 
in the model, and there were no piezometers downstream from leaf. 

All pressure measurements were made with water-filled U-tubes, and the 
readings were made in tenths of an inch and estimated to the nearest hun- 
dredth of an inch. ‘ The test procedure for each gate consisted of setting the 
leaf to the desired position, allowing a few minutes of operation for conditions 
to stabilize, and then taking the pressure readings. The leaf was then set to 
the next desired position and the procedure was repeated. 


Investigations 


The downpull forces on the gate leaves were determined by measuring the 
pressures acting on the top and bottom surfaces of the leaves, and not by 
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measurement of stem loads by strain gages, spring scales, weights, or other 
devices. Previous experience at the Bureau of Reclamation has shown that 
the forces determined by pressure distribution may be more accurate than 
forces determined by stem load measurements because of the effects of fric- 
tional factors, vibration, and leaf movement. Proper clearances between 
parts are of particular importance, and it is often difficult to predict what 
these clearances will be, or to produce them in the model. When using the 
pressure-area method, care must be taken to insure getting pressure read- 
ings in all regions which can effect the downpull on the leaf, and to apply 
these pressures properly. 

In the case of the Glendo fixed-wheel gate, the only condition where down- 
pull will be an important factor is during submerged operation when gate 
Openings are small relative to the leaf thickness. In this case the only pres- 
Sure relief obtained beneath the leaf is that due to circulation of water be- 
tween the leaf and the jet. 

The pressures on top of the leaf can vary from essentially atmospheric 
when the backwater is insufficient to fill the tunnel, to high values when the 
downstream tunnel is flowing full. 

When no backwater is present and the gate opening is appreciable, the jet 
will pass beneath the leaf and clear the downstream web enough to allow free 
aeration. Thus, the pressures under the leaf will be about atmospheric. The 
pressures above the leaf will also be about atmospheric because the gate 
seals are on the upstream face and the bottom surface is open to the tunnel 
pressure at the downstream face (Figure 1). There is little or no downpull 
under these conditions. 

In the slide gate, downpull will be an important factor through a wide range 
of openings for both free discharge and submerged operation. The maximum 
downpull can occur with submerged operation because subatmospheric pres- 
‘sures can exist on the flat-bottomed portion of the leaf during this operation. 
In contrast to this, the lowest pressures on the flat-bottomed portion of the 
leaf during free discharge operation are about atmospheric. The bonnet pres- 
sures will be high at small gate openings for either free discharge or sub- 
merged operation because this type of gate seals on the downstream face of 
the leaf, and is open to the tunnel pressures at the upstream face (Figure 2). 
‘The high bonnet pressures are the main factor in the high downpull forces 


produced on these gates. 
¢ 


-Fixed-Wheel Gate 

In the simplified air model equivalent values of head differential that will 
occur across the full-size gate could not be set. But, as in a water model, 
the pressures on and across the leaf varied in accordance with the model 
heads and velocities. Thus, prototype pressures were obtained by multiplying 
‘model pressures by the ratio of the computed prototype head differential 
across the gate to the measured model differential. To compute the prototype 
head drop, it was assumed that the 21-foot-diameter tunnel was ruptured at 
e power bifurcation, and that the reservoir was at the maximum water sur- 
face elevation of 4669.0, producing a head of 178 feet. The head loss through 
e intake structure and tunnel with the gate 100 percent open limited the 
aximum flow to 21,500 cfs. As the gate closes to reduce the flow, the tun- 
1el losses become smaller and the head drop across the gate becomes larger. 
very small openings nearly the full shut-off head of 178 feet is acting on 
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the gate. The relation of prototype head drop to gate opening is shown in 
Figure 6A. 

The model was tested at a number of gate openings, and the pressure drop 
across the gate and the pressures on the bottom of the leaf were measured at 
each of the openings. These measured pressures are shown in Table 1. It is 
interesting to note that the pressure just above the lower seal (Piezometer 
No. 1) is almost the same as the pressure acting on the roof at the leaf, and 
hence on the upper seal. The upper seal has the same plan area as the lower 
one, and therefore the upward force on the upper seal is the same as the 
downward force on the lower seal. These forces balance one another; thus 
Piezometer No. 1 is not used in the downpull calculation. From these data, 
the prototype downpull forces were calculated as follows, using as an example 
a 3-inch prototype gate opening. 


Model pressures in inches of water were: 


Upstream Downstream Gate web ; 
head head eat (average) pitas 
+8.07 -0.09 -1.25 -0.51 -0.52 

Prototype head drop (from Figure 6A) = 178 feet 
; 178 
Factor for conversion = 3.07 + 0.09 = 21.80 


Downstream head = -(0.09 (21.80) = -1.96 = head on top of leaf 
Under the seal = -(1.25) (21.80) = -27.21 

Under the web = -(0.51) (21.80) = -11.12 

Under the DS lip = -(0.52) (21.80) = -11.34 


If it is assumed that the pressure distribution at any section along the leaf 
is the same as the measured distribution at the center line section, the down- 
pull may be found by multiplying the pressure difference between the gate top 
and bottom by the specific weight of water and by the area of the gate section 
(Figure 6B) on which the pressure difference acts. 


Table 1 


MODEL PRESSURES(IN INCHES OF WATER) ON GLENDO FIXED-WHEEL GATE 
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Seal = (27.21 - 1.96) (62.4) (16.5) (0.3802) = 9,900 pounds 
Web = (11.12 - 1.96) (62.4) (16.5) (2.9282) = 27,610 
DS lip = (11.34 - 1.96) (62.4) (16.5) (0.2240) = 2,160 

Total 39,670 pounds 


A plot was made of the computed prototype downpull forces acting on the 
seal, under the lower web, and under the downstream lip for a wide range of 
gate openings (Figure 6B). The total downpull, which is the sum of the 
individual downpulls, is also shown. 

This downpull information may be applied to other geometrically similar 
installations by the following method: 


A. Find the downpull that would occur at the desired gate openings with the 
head drop for the proposed installation: 


Hy 
DP; Drie 


DP; = Downpull, Figure 6B 
DP» = Downpull with H, 


jen) 
ue) 
Ul 


Head drop across Glendo gate, Figure 6A 
H, = Head drop for proposed installation 


B. Obtain the desired downpull by multiplying DP» by the ratios of the 
areas under the two gates: 


D An 
P3 = DPy x i 


DP3 = downpull under the proposed gate with Hy 
Ag = area under the tested gate 
A, = area under the proposed gate 


It will be noted in the sample computation that the pressure beneath the 
bottom seal reached a dangerously low value of -27.21 feet of water. This is 
low enough to produce cavitation in installations at fairly high elevations. 
Even at elevations as low as sea level, pressure fluctuations of only moder - 
ate proportions would be enough to lower the pressures momentarily to where 
cavitation would occur. If the 178 foot head differential across the gate were 
increased to 200 feet, the seal pressure would become -30.60 feet and cavita- 
tion could extend over much of the seal. At a 250-foot head differential, cavi- 
tation pressures would exist for gate openings from about 2 to about 4 inches. 
For the conditions where cavitation is general, downpull on the seal is com- 
puted by applying cavitation pressures to the area under the seal. 


Slide Gate 


It has been previously noted that a model that uses air as the flowing fluid 
and that discharges into the atmosphere is operating under submerged 
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conditions. In the case of the slide gate, the flow patterns along the upstream 
face of the leaf and along the sloping bottom of the leaf will be the same for 
free discharge and for submerged operation. Pressures measured on these 
components in an air model are therefore equally applicable to free discharge 
and submerged conditions. A distinct flow difference will exist, however, 
along the horizontal or flat portion of the leaf which lies downstream from the 
abrupt ending of the slope. During free discharge the flow will spring clear 
of this surface and the pressure on it, and on any seal extension, will be about 
atmospheric. During submerged operation an eddy of the flowing fluid will be 
in contact with the surface and the pressures will be a function of the head, 
the gate opening, and the effective back pressure on the gate. The lowest 
pressure that can occur on the surface when water is flowing is the vapor 
pressure of water, about -30 feet gate. The pressure under a seal extension, 
if one is used, will be affected in about the same manner. 

Practical considerations dictate the width of the flat on the leaf bottom and 
on any seal extensions required. Therefore the geometry of the gate bottom 
may not be exactly similar for gates designed for different sized installations 
and heads. To make the data applicable to these variations in design the data 
was reduced to that applicable to the sloping portion of the leaf, and that ap- 
plicable to the flat-bottom portion, including the slot projections and seal ex- 
tensions (Figure 84). 

The method of determining the slide gate downpull was as follows: The 


‘model bonnet pressures, and hence the pressures acting downward on top of 
_the leaf, were divided by the head differentials across the gate. This dimen- 


sionless ratio was plotted against percent gate opening in Figure 8B. Simi- 
larly, the pressures acting upward on the sloped portion of the leaf bottom 
were divided by the total heads producing flow. This ratio was also plotted in 
Figure 8B. The difference between these ratios, at any gate opening, multi- 
plied by the prototype head producing flow, equals the prototype differential 
head acting downward on the cross-sectional area of the sloped portion of the 
leaf. This head difference, multiplied by the appropriate area, and by the 
density of water, results in the downpull force on this part of the leaf. To 
this force, the forces acting on the flat portion of the leaf and on the slot 
projections and seal extensions are added. This sum equals the total down- 
pull force. 

The pressures acting on the sloping portion of the leaf were investigated in 
detail. Three rows or piezometers were included on the leaf, one 0.18 inch 
from the leaf side, the second 1.10 inches from the leaf side, and the third in 
the leaf center line (Figure 5). The pressures for each row are plotted non- 
dimensionally in Figure 7 for gate openings of 10 through 80 percent. The 
pressure at any point on a similarly shaped leaf bottom can be determined 
by multiplying the appropriate dimensionless factor by the head producing 
flow (H~-hg) and subtracting this result from the total head upstream from 
the gate. For determining the downpull forces, the model pressures were 
replotted dimensionally with the forces being considered vertical. Planimeter 
measurements were made of the areas within the pressure distribution en- 
velopes for the sloping portions of the leaf. These areas were divided by the 


2.23 inch thickness of the leaf above the sloped bottom to give an average 


pressure value at each piezometer row. The average pressures obtained for 
the row 0.18 inch from the leaf side were assumed to act from the leaf sides 
to stations 0.50 inch from the sides. The pressures for the row 1.10 inches 
from the side were assumed to act from the 0.50 inch stations to stations 2 


January, 1959 


HY 1 


04 


V3GOW YIV WONs viva 


4VI1 31V9 30INS 40 WOLLOG NO SAYNSS3IYd-2 JYNDIS 


YIU} 8b X @PIM,1Z'@ Oa) japoy 


SSOUDIY} $09] =) 


309| JO 890} Woad4sdn Woy aouD\sIG =x 


SNIN3dO 3Lvd %08 


o> 4v3a0 30/S 3¥37 WOU Or 
or 80 90% 75.5 Mn 
i Pe le, Teel 
b 
b= 4 
= 
a | oaks 
QNIN3dO 31V9 %09 
¥4v37 3018 wan WOUd Or 
0190 90°%»n 20 0 
if 
b- 


ONIN3dO 31LV9 %Ob 
dvavnI 30'S 4¥37 WOU ,01'I 


---uo}pouado ebsoyosip 

884) YIM UOI4Jod 

WO440q-40)) UO aunssoig 

Ll 

ONINAdO 31lvd %02 
30ls aay WOUd 01") 

Sy tJ 


301S 4V37 WOU 810 
O1_ 80. 90 "po 20 


3G1S 3¥31 WOUS ,8I'0 
90 a 


0 900 20 


30'S 4¥31 WOHS ,81'0 
als 


aoceeaemersers 0140480 


U WO}40q-4D|J UO ass aig 


31S 4¥37 WOH 810 


AM 
90 


4H 


@y-4y 
Xy- 


=f 


Zy-4y 
Xy-1y 


[pabsowgns yim uo}4s0d 


Zy-1y 
Xy-4yq 


3406 Wo., woasysumop poay asnssaig = 2y 


S48}9WOZAId joa] }o POaY. auNssauy = XY 
9406 wosy wWoasysdn poay jojo) = ty 


SNINAdO 31V9 %O2 
>4vaq 30'S 4V37 WOU 01"! 


ol! 80 


30'S 4V37 WOUS , 


ol 90 od %y0 20 


ONIN3dO Jlv9 %OE 


dsav70 30'S 4¥37 WOU, OI 
, din oe ah, 
0. a2 
b 
f~ 
SNIN3dO 3LV9 %Ol 
3437 


30'S 4¥37 WOUd,,01"I 
gore 


0 


45-409) 40 Uo1Juod 
i PaWO440q 4014 


4-409) 40 UOI,J0d Buidoys- 


Se ee ee ee a ee 


31S 4V37 WOYS ,8!'0 
alle, 


Xy-4y 


30!S JV37 WOUS,,8I'0 


oy 


30'S 4¥37 WOUd ,BI0 


al/&, 
0 


if 
pix 
30'S 4¥31 WOUS ,81'0 
mn 2 
a |) 
409) 40 zr|=x 
9904 woes4sdr)-.._ ; +4 
wl 


coeeneeneer== (D8) 40 
990) woa.4sumog 


ASCE HYDRAULIC DOWNPULL 55 


inches from the sides. The pressures at the center line piezometers were 
assumed to act over the remaining width of the leaf. The products of these 
average pressures and appropriate areas were added together for each gate 
opening. These pressure-area summations were divided by the projected 
area of the sloping portion of the leaf to give the equivalent single pressure 
values that can be assumed to act over the sloped area of the leaf at each 
gate opening. These pressure values, divided by the total head producing 
flow, establish the curve plotted in Figure 8B. At zero opening, the average 
pressure acting upward on the plan area was the same as the total head. It 
decreased to a minimum at 70 percent open and increased again at 80 percent 
open. No tests were made at openings above 80 percent because the downpull 
would be small. 


Free discharge conditions 

In the full-sized slide gates the back of the leaf carries a seal that is in 
contact at all gate openings with a seal plate in the bonnet. There can be no 
flow or appreciable leakage between the back of the leaf and the bonnet, and 
only a little flow down the gate slots. Thus the bonnet pressure is about the 
Same as the stagnation pressure just ahead of the leaf. The model pressure 
measurements, based on the piezometer in the top of the conduit 1 inch 
ahead of the leaf showed that without submergence the bonnet pressure was 
the same as the total head in the conduit at zero gate opening, and that it de- 
creased with respect to the total head as the gate was opened (Figure 8B). 

The difference between the bonnet pressure and the average pressure Over 
the sloped bottom was a maximum at an opening of about 45 percent and the 
difference did not exceed about 53 percent of the total head (Figures 8B and 
C). The net downpull on the flat bottom, the slot projections, and any seal 
extensions can be determined from their areas and the difference between 
the bonnet pressures and atmospheric pressure. The total downpull will be 
the sum of the individual downpull forces acting on the various portions of the 
leaf. 

A sample calculation of the downpull forces follows: 


Assume: Free discharge conditions, gate leaf 6.0 feet wide, 2.0 feet thick, 
3-inch flat on bottom, 6- by 6-inch slot projections, no seal ex- 
tension (Figure 8D) Head H = 200 feet. 


The net force acting downward on the plan area over the sloped part of the 
leaf bottom is: 


TO 
FES | 


: Tf 
2) | Head di To 
HT 


2) x HT = head Pie) 94,2 105.0 
diff (feet, H2O) 


pete aati 


3) x 62.4 = 2,658 | 5, 878 aren 

Press diff (psf) 
d 4) x area = 27,909 | 61, 719 | 68 69, 321} 6 : 
‘ (5) | downpull, 1b 
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The net force acting downward on the plan area of the flat bottom and the 
slot projections (no seal extension on standard gate) is: 


Bonnet Pres %] 99.5 06.0) Sao lips DOE 
184, 172 


HT 
(6) x HT = 199.0] 192.0 
Bonnet Pres = 
Head diff (ft, 
H20) 
Gi) x 62.4, = 12, 417 | 11, 981 (ae 507 | 10, 770 


(8 
Pres diff (psf) 

~(9)] (8) x Area = 24, 834 | 23, 962 | 23,014 | 21, 540 
Downpull; A = 
(0.25 x 6.0) + 
20.5 x 0.5)5= 
2.00 

(5) + (9) = 52, 743 | 85, 681 | 91, 810 | 90, 861 
Total downpull 
lb 


(i0 


(Hy = Total head upstream of gate.) 


If it were assumed that the full total head acted over the entire sectional 
area of the leaf, the calculated downpull would be 200 x 62.4 x (10.5 + 2.0) 
= 156,000 pounds. The maximum downpull of 91,810 pounds obtained with the 
above test information is 59 percent of this value. This percentage will vary 
somewhat depending upon the width of the flat bottom of the leaf, the location 
of the leaf projections, and whether or not seal extensions are used. 


Submerged Conditions 

If the tailwater elevation rises to submerge the gate or fill the downstream 
conduit, and the headwater elevation remains unchanged, the head differential 
across the gate, and hence the rate of flow, decreases. In spite of these 
changes, and provided that the Reynold’s Number remains reasonably high, 
the flow pattern approaching and passing along the face and sloping portion 
of the leaf remains the same. And because the flow pattern is the same, the 
difference between the upstream conduit total head and the head at any point 
on these areas of the leaf, divided by the head differential across the gate, 
peeins the same for any gate opening. Expressed mathematically, 
{a5 - i = K, and the downpull force on the sloped portion of the leaf will be 
the same as for the free discharge condition with the same differential head. 
In the case where a gate can be placed more deeply or less deeply below the 
tailwater elevation, the effect of submergence change is merely to equally 
raise or lower the bonnet and leaf slope pressures. The head differential 
across the gate would remain constant, and so would the downpull forces 
(assuming no cavitation at the lower submergences). 
When the discharge is into a pool, the head producing flow may be taken as 
the difference between the total head upstream from the gate and the pool 
depth downstream. When the discharge is into a filled conduit, the head 
producing flow may be taken as the difference between the total head upstream 
and the static head a distance of about 3 conduit heights downstream. 
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For submerged flow where conditions would permit vapor pressure under 
the flat on the leaf, the method of calculation would be the same as that in the 
free discharge example except that the downpull head differential in Item 7 
would become the total bonnet pressure plus the numerical value of the vapor 
pressure, about 30 feet of water. This produces the greatest downpull pos- 
sible, but it is unlikely that vapor pressure would ever occur over this whole 
area. Thus, this extreme downpull condition will probably never be reached. 
In the more likely cases where conduit failures or unexpectedly low tailwater 
will not be factors, and where the submergence will always be sufficient to 
hold the pressures on the leaf flat and on the seal above vapor pressure, the 
net downpull would be determined from the difference between the bonnet 
pressure and these surface pressures. Such surface pressures may be diffi- 
cult to predict exactly, and would probably have to be determined experimen- 
tally if precise results were needed. But for ordinary work, the pressures 
can be obtained from the non-dimensional data presented in Figure 7. 


SUMMARY 


Air models are reliable tools for determining flow phenomena and pres- 
sure distribution in fluid flow passages. As such, they are readily adaptable 
to studies of downpull forces on gates, particularly where the gates operate 
submerged. They may also be used for downpull studies when the gates dis- 
charge freely, provided that care is taken in analyzing the flow conditions and 
resulting pressure distributions on the various portions of the gate leaves. 
Air models are less expensive to construct and test than hydraulic models, 
and they are accurate, convenient, and easy to work with. 

The two gate types discussed in this paper are typical of installations bein: 
made throughout the country. The data are presented in a form usable for 
any sized gate, provided that the top and bottom of the leaves are generally 
geometrically similar to the test gates. The maximum downpull on the tested 
fixed-wheel gate occurred during submerged operation at about a 3-inch open 
ing. The downpull rapidly decreased as the opening increased. On the slide 
gate, the maximum downpull force occurred at about a 45 percent gate open- 
ing. This calculation assumed that the total head remained constant at the 
gate regardless of gate opening. In the usual case, friction in the conduit 
ahead of the gate will cause a loss of head as the flow rate increases. The 
maximum downpull would be less than in the example, and would occur at a 
smaller gate opening. Bonnet pressures play a particularly important part in 
the case of the slide gate, and in designs where leakage from the bonnet is 
permitted or encouraged to decrease the bonnet pressures, downpull will be 
materially reduced. Of course, this drainage or leakage may create other 
complications, and cannot be regarded as an automatic solution to the down- 
pull problem. 
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APPLICATION OF SNOW HYDROLOGY TO THE COLUMBIA BASIN 


Oliver A. Johnson, 1 M. ASCE and Peter B. Boyer2 


SYNOPSIS 


Observations of the weather elements that cause the melting of a snow 
pack were made for several years in three laboratory areas located in dif- 
ferent weather regimes of the high mountains of the western United States, 
and the data obtained were analyzed and correlated. Starting with the inter- 
relationships between the weather elements that produce snowmelt at a point, 
methods and procedures were found for the rational analysis of snowmelt in a 
relatively small basin. The application will be enlarged to cover large parts 
of the Columbia Basin as the network of stations covering the basin is en- 
arged and additional climatological phenomena are measured. 


INTRODUCTION 


The Columbia River is principally a snow-fed stream. Its average annual 
unoff ranks very high among the streams of the United States, averaging 
bout 180,000,000 acre-feet. The basin is located just off the Pacific Ocean 
n the belt of the heavily moisture-laden prevailing westerly winds. In winter 
he Aleutian lows attract maritime air masses which deposit their moisture 
iS they ascend the high mountain slopes on the west side of the Cascades. The 
astern Slopes of the Cascades and the plains of the Columbia Plateau lie with- 
n the “rain shadow” of the Cascades, and annual precipitation in this region 
nay be as little as eight inches. Additional uplift of the air masses by in- 
erior mountain ranges and the Rocky Mountains causes increased precipita- 
ion in the headwater regions, in amounts ranging to more than 80 inches per 
ear. 
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The Columbia rises in Columbia Lake at elevation 2652 feet m.s.l. as 
compared to elevation 1670 for Lake Itasca, the source of the Mississippi. 
The mean elevation of the Columbia drainage east of the Cascades is about 
4500 feet m.s.l. The rugged topography is important because of the large 
vertical distance through which the water descends on its way back to the 
Pacific, and because of the possibilities for economical construction of large 
dams afforded thereby. In January 1957, 25 million kilowatts of hydro capaci 
ty were available in the United States while that in Oregon, Washington, Idaho 
and that part of Montana west of the Continental Divide (comprising about 9 
percent of the U. S. area) totalled more than 7 million kilowatts or 28 percent 
of the total. 

During the design of a storage project on the Columbia or one of its tribu- 
taries, snow hydrology must be considered in the computations required in th 
selection of the project design flood. This is the largest flood that must be 
controlled if a satisfactory flood control plan is to be achieved. For that par 
of the Columbia and its tributaries lying to the east of the Cascade Range, the 
melting of the snow pack creates the severest design conditions. In the Willa 
mette and other coastal streams the most intense floods are caused by sever 
rain-on-snow events occurring in the period November through February. 

Power projects and flood control storage are entirely compatible in the 
storage projects above Bonneville Dam. Power draft through the fall and 
early winter will provide ample space for storage of floodwaters from moder 
ate rainstorms that have occurred infrequently in December, and continued 
draft will provide additional space before the onset of the snowmelt season. 

In the storage projects above Bonneville Dam there are monetary benefits 
to be derived by releasing stored water for power production at site and at 
downstream plants, and there are monetary benefits derived from flood 
damage prevented by such storage, both locally and in the vulnerable reach 
of the Columbia lying downstream from Bonneville Dam. 

In the operation of the Columbia River projects, the principal objective is 
to provide storage space for a possible snowmelt flood and still have full 
reservoirs at the beginning of the power draft period. To achieve this ob- 
jective, the water content of the snow on the ground is measured at frequent 
intervals throughout the accumulation and depletion period. With the current 
observations of precipitation and temperature and assuming normal weather 
for the remaining period, satisfactory runoff volume forecasts are made. 
These forecasts are increasingly more accurate as the season progresses 
because the future precipitation becomes progressively less important. Larg 
deviations from normal precipitation during the melt period in May and June 
can cause the runoff volume to differ considerably from the forecasted event 
but operating parameters can be devised that will afford satisfactory regula- 
tion of the projects. 

There are other snow-fed streams that are of great importance to the 
United States, such as the Colorado, Sacramento, and Rio Grande, to name a 
few on which the Corps of Engineers and Bureau of Reclamation have built 
many reservoirs. A knowledge of the rates of accumulation and depletion of. 
snow are of utmost importance for intelligent and efficient operation of these 
reservoirs. The Weather Bureau, Soil Conservation Service and Geological 
Survey also have a great interest because they are the principal agencies 
charged with the collection, processing and dissemination of hydrologic in-- 
formation. All are endeavoring to improve their forecast of the available 
basin snow water, of its rate of release and of losses. Reliable forecasts by 
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elatively rapid means is mandatory for both economical design and opera- 
on. 

The need for detailed, accurate knowledge of the characteristic and be- 

avior of snow was recognized by several of the Federal agencies and after 
me discussion, a program known as the “Cooperative Snow Investigations” 
aS agreed upon and supported by the Weather Bureau and the Corps of Engi- 
eers in 1945. Limited assistance was furnished by the Geological Survey, 
reau of Reclamation and Forest Service. Direct participation by the 
Weather Bureau ended in 1952, and thereafter the basic research was con- 
tinued by the Corps of Engineers under the name of “Snow Investigations.” 
Throughout its life the program was pointed toward the solution of hydrologic 
problems pertinent to mountain regions of western United States. Three 
laboratory areas were heavily instrumented and studied intensively for peri- 
ods of five to eight years in order to provide the basic data for determination 
of the hydrologic factors involved. These laboratory areas were designated 
Central Sierra (near Soda Springs, California), Upper Columbia (near Marias 
Pass, Montana), and Willamette Basin (on Blue River, in Oregon). 

Beginning in 1954 a part of the group effort on snow investigations was 
directed toward compiling the results achieved by these studies in the form 
of a Summary report that would serve as a reference volume for hydrologists. 
The work was concluded in June 1956 with a volume titled “Snow Hydrology” (1) 
containing some 400 pages of text and 72 plates. The principal features that 
apply to the Columbia River Basin are presented in this paper. 


Point and Areal Relationships 


The field of snow hydrology is concerned with the evaluation of various 
components of the hydrologic cycle. It draws upon the sciences of meteorolo- 
sy, hydraulics, thermodynamics, geology, soil mechanics and plant physiology. 
The development of the applied snow hydrology begins with observations and 
measurements, at a point, of factors affecting each component. This is fol- 
owed by an analysis, leading to the establishment of mathematical relation- 
ships. Prior to basin wide application of the point values, the relative magni- 
udes of the components are determined, and emphasis is placed on the more 
mportant elements. Having determined point values of the components of the 
tydrologic cycle, the next problem is to modify the equations to fit the local 
slimate and topography peculiar to an individual project basin. During the 
‘ourse of the snow investigations it was discovered that experimental data 
rom areas as large as the laboratories (4 to 20 sq. mi.) could be fitted to 
oint relationships only with difficulty, and near the end of the program ex- 
erimental data from impervious areas measured in the hundreds of square 
eet were obtained to assist in evaluating the coefficients needed to determine 
int melt. 

Having determined point values, the next step consists of utilizing these 
alues to determine amounts and distribution over a basin area, with respect 
0 differences in environment, as well as to time. This involves procedures 
yhich are much less exact than those for point evaluations. It is impractical 
0 attempt a point-by-point analysis in basin application; rather, it is neces- 
ary to deal with basin averages in major subdivisions of geography or en- 
ironment, and also to deal with averages in time. This concept leads to the 
se of point measurements as indexes to represent basin values. Since the 
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use of indexes involves the theory of sampling errors in measurement, statis 
tical procedures may be employed to evaluate the reliability of estimates and 
also the weights assigned to individual factors. Intelligent selection of indexe 
requires an understanding of the heat processes and the relative effects of 
topography, forest and sky cover on the heat supply to the snow pack at a 
point. 


Water -Balance in Areas of Snow Accumulation 


When reliable hydrologic data are available, but the records are too short 
to justify a statistical approach, an accounting method, called water-balance, 
is employed in forecasting the stream runoff during a specified period. In 
addition, a knowledge of the water-balance for a given area is necessary in 
order to select appropriate forecasting parameters and to interpret their re- 
liability. The water-balance technique may be used as a forecasting proce- 
dure by quantatively balancing runoff against precipitation, change in snow- 
pack water equivalent, and losses. A knowledge of the water balance is also 
necessary for rate-of-flow forecasting and for design flood computations. 

The water-balance method is formulated by the equation: 


Qe Wye Pew, an (1) 


in which P is the net basin precipitation (after correcting for gage catch de- 
ficiency and interception loss); W and W9 are the basin snowpack water 
equivalents at the beginning and end of the period; and L is the loss con- 
sisting of evapotranspiration (Let) and soil moisture change (Lg), that is, 


L = Let # Ls (2) 


in which Lg is positive when water is added to the soil to reduce its moistur 
deficiency and negative when moisture is being evaporated from the soil; Q 
is the generated runoff volume during the period considered, equal to the ob- 
served runoff corrected for recession flow. Generally the purpose of a run- 
off forecast is to determine in advance the volume of water (Q) expected to 
enter the reservoir between two given dates. The forecasted Qg is therefore 
modified to reflect this desired volume as follows: 


Ore Ot 0, oa i (3) 
whence equation (1) reduces to 


QaW, 4 P-Wo-L4AQ-% (4) 


in which Q; and Q9 are the runoff volumes under the initial and terminal 
recession curves for the period considered. 


Snowmelt Relationships 


Heat is required to melt the snow. At present we know enough about the 
principles of heat exchange between the snow and its environment to formula 
general snowmelt equations based upon detailed studies for each of the three 
laboratory areas. Furthermore, these equations can be expressed in terms: 
of the ordinarily available meteorological elements which properly represen 
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the heat processes. The level of our confidence is not as high in the applica- 
tion of the point snowmelt principles to basin-wide melt, since only a small 
number of basin-wide melt analyses have been completed. Each basin pre- 
sents a separate snowmelt problem and has a unique solution. The principal 
problem is to integrate the effect of topography on the climate near the snow. 
This involves (1) the separation of the project basin into areas of homogeneous 
| heat supply under cloudless, cloudy, and rainy sky conditions, and (2) deter- 
mination of the basin constants for each area—constants which will integrate 
the influence of topography on solar radiation, on long wave radiation, and on 
convective-condensation heat exchange. 
The processes responsible for the supply or less of heat are: 


Absorbed solar radiation (H,g) 


Net longwave radiation exchange between the snowpack and its environ- 
ment (Hy]) 


Convective heat transfer (sensible heat) from the air (He) 
Latent heat of vaporization released by condensate (He) 
Conduction of heat from underlying ground (Hg) 

Heat content of rain water (Hp) 


These processes are influenced by meteorologic, geographic and time factors 
which must be considered and related to the ordinarily observed data. The 
general equation for snowmelt at a point is: 


M:- ZH = £H = 0.0051£H (5) 
203B 203(0.97) 


in which M is in inches, 2H is the algebraic sum of heat energy by short 
wave radiation, long wave radiation, convection, condensation, conduction and 
precipitation; 203 (= 80 x 2.54) is the number of calories required to melt one 
inch of water equivalent of ice at 0° C over a unit area; B is the thermal 
quality which is equal to the ratio of heat required to melt a unit weight of the 
snow to that of ice. B>1 when the snow temperature is less than 0°C and 

B = 1 when the snow is at 0°C and dry. But a melting snow, after drainage, 
contains liquid water which requires no heat to melt. Therefore, its thermal 
quality, after drainage, varies somewhat with density. An average value of 
B= 0.97 is satisfactory, as shown in equation (5). 

The analysis of observations and measurements of radiation, temperature, 

humidity, condensation (or evaporation) and melt at a point resulted in the 
following equations, where snowmelt, M, for each heat component is in inches 
per day. 


Shortwave Radiation Melt 


Le = 0.0051 (l-a) I; (6) 


where a is the albedo (reflectivity) of the snow surface, expressed as a deci- 
‘mal fraction and 


I, is the incident shortwave radiation in langleys (Cal. / sq. cm.). 
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Longwave Radiation Melt 
For clear skies 


M., = -022 (T, - 32) - 0.84 (7) 


where T, is the air temperature at the 10-foot level in OF, 
For cloudy skies (low clouds) 
M_, = 0-029 (T, - 32) (8) 
where T,, is the temperature of the cloud base in aw 
For forested area 


Ms 0.029 (T, - 32) (9) 


Convection Melt (theoretical) 
-1/6 
M, = 0.00629 (p/p,) (Z,2,) (rey. (10) 
where p/ Po is the ratio of the atmospheric pressure at sk site to standard 
sea-level pressure. 


Z, is the height in feet above the snow surface of the air temperature 
measurement. 


Zp is the height in feet above the snow surface of the wind speed measure- 
ment. 


T, is the air temperature in OF. 

T, is the snow surface temperature. 

Vp is the wind speed in miles per hour. 
Condensation Melt (theoretical) 

M, = 0.054 (Z, 2)"1/®(e, - .) v, (11) 
where Z, is the height in feet above the snow surface of the vapor pressure 
measurement. 

Zp and Vp are as defined above. 
€g is the vapor pressure of the air in millibars. 


e, is the saturated vapor pressure of the snow surface. 
Combined Convection and Condensation Melt (simplified) 
= / 0.0084 v 7 [0.22(T, - 32) # 0.78 (Tq - 32)7 - “(237 
where v is the wind speed at the 50-foot level in miles per hour. 


T, is the air temperature at the 10-foot level in OF. 
Tg is the dew point temperature at the 10-foot level in °F. 
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Conduction of Heat from Ground 


Melt from this source is negligbly small, so that direct evaluation is un- 
necessary. Average value is approximately 0.02 inches/day. 


Rain Melt 
M = 0.007 (T.. - 32) P.. (13) 


where 7. is the temperature of rainwater in °F (usually assumed to be that 
of the air temperature at the 10-foot level). 


>. is the rainfall in inches. 


A practical method of establishing a reliable clear-weather melt equation, 
in terms of the available heat indexes, for a ripe snow in a given basin is by 
means of statistics. The influence of topography and forest cover on each 
heat process will then be integrated in derived coefficients that will also in- 
clude the period of melt and the constants of proportionality between units of 
measure. The proposed form of such a linear regression equation is: 

Q, = k'(1-a)I, # AT, k,, (DT, # CT, -T,)v 4 D (14) 
where Qg is the generated runoff. The term, k! (1-a)l;, is the solar radia- 
tion melt in which Jj is the solar radiation incident on the horizontal surface 
of the snow area, a is the albedo of the snow, and k' is the basin constant 
integrating the effect of topography and forest cover on solar radiation. The 
term CT, defines the long wave radiation melt with C as its derived coef- 
ficient. The term, kee(BT,+ CTg - Tp) v, is the convection-condensation 
melt, in which Keg is a constant defining the degree of basin exposure to 
wind; T, and Tg are the air and dew point temperatures; B, C, and Tp are 
derived constants, and v is the wind speed. D is the regression constant. 
Ground melt is negligibly small and it will be partly reflected in the other 
terms. 

The melt at the snow surface for the generation period is: 


Wo. 07 1 (15) 
g 

in which L is the loss of water by evapotranspiration and deep percolation. 

Note that the proposed statistical analysis assumes L as a function of T, Tg 

and I. Evapotranspiration is a function of these elements; and deep percola- 

tion, unless definitely established, is assumed to be negligible. 

Examination of data taken on cloudy and cloudless days indicates that an 
adjustment factor for average cloudiness should be used in computing long 
wave radiation melt. On the days of complete cloud cover and on rainy days, 
it is easier to develop surface melt equations by a rational approach. 

Melt equations peculiar to each sub-basin may be derived statistically by 
direct correlation of appropriate meteorologic factors to daily generated run- 
off amounts as explained above. In many cases, it may not be feasible to de- 
rive basin snowmelt equations in this manner because of lack of appropriate 
data. General snowmelt equations were therefore developed from laboratory 
data, which may be used to evaluate basin snowmelt on the basis of the local 
physical characteristics and the meteorologic factors affecting snowmelt. The 
physical characteristics of a basin affecting heat supply may be expressed in 
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terms of coefficients of average forest cover, exposure to wind, and exposure 
to solar radiation. The form of the equation is determined by the amount of 
forest cover, since it has a marked effect upon the relative importance of 
both wind and solar radiation in estimating snowmelt. For rain-free periods, 
the meteorological parameters include free air and dewpoint temperature, 
wind speed, solar radiation and cloud cover. During periods of active melt, 
the only variable of the snowpack condition is the albedo and the snow cover. 
For estimating snowmelt during rain, values of air temperature, wind speed, 
and rainfall are required. 

The components of the clear weather melt (in inches per day) were each 
modified by the appropriate basin factor to account for slope and orientation 
(k'), forest cover shading (F), exposure to wind (k), and for cloud cover (N), 
resulting in the following equations: 


Shortwave radiation melt, Mys: 


M_, = 0-0051 k'(1-F)(1-a) I, (16) 
Longwave radiation melt, My}: 
M y= (1-N)(1-F) (0.022 (1, -32) - 0.84 74 (N # F(1-N}7 0.029(T,-32) (17) 
Convection-Condensation melt, Mee: 


Mog = 0-0084 kV5, (0.22 (T, -32) ¢ 0.78(Tg -32)_7 (18) 


Ground melt: a constant 0.02 in. per day. 


In the above equations Ij is the daily short wave radiation (in langleys) 
reaching a horizontal surface at the latitude and altitude of the snow area; 
Ta and Tg are the daily average air temperature and humidity at 10 feet 
above the snow; ana V5 is the daily average wind speed (in mph) 50 feet 
above the snow. This study was followed by a statistical analysis of the clear 
weather snowmelt in each of the three laboratories, using estimated basin 
average temperature, wind, k, and F values. The analysis resulted in the 
following typical snowmelt formulas, in which k, k! and F are estimated from 
a careful examination of the topography of the snow area. 

Heavily forested area (like Willamette Basin Snow Laboratory) having 
more than 80 percent forest cover: 


M = 0.074 (0.537, 40.47 T, ) (19) 


Forested area (like Skyland Creek, Upper Columbia Snow Laboratory) hav- 
ing 60 to 80 percent forest cover: 


M=k (0.0084 v)(0.22T, # 0.78 Ty) 4 0.0297, (20) 


Partly forested area (like Central Sierra Snow Laboratory) having 25 to 
60 percent forest cover: ; 


M = k'(1-F) (0.00401; ) (1-a) # k(0.0084v) (0.22T! ?# 0.78T3) ? 


21 
F(0.029T') Me 


Open area (less than 25 percent forest cover): 
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M = k'(0.00511;)(1-a) # (1-N)(0.022 TL - 0.84) 4 N(O.029T!) ¥ 


k(0.0084 v)(0.22T! ¢ 0.7874) 
here: 
M is the snowmelt in inches per day over the snow covered area. 


t4 is the difference in air temperature measured at 10 feet above the 
snow and the melting snow surface temperature, in °F. The fore- 
casted temperature of the air at a station may be extrapolated to the 
pow level by the lapse rate. The temperature of the melting snow is 
32 °F. 


T is the difference between the dewpoint temperature at 10 feet above 
the snow and the melting snow surface, in °F. It is approximated 
from upper air or surface humidity forecast. The dewpoint of the 
melting snow is 32°F. 


f is the difference between the cloud base and the melting snow surface 
temperatures in °F. It is estimated from the upper air temperature 
or by a lapse rate from the surface temperature. 


v is an estimated representative wind speed at-50 feet above the snow, in 
miles per hour. Wind stations in the proximity may serve as a guide. 


is the estimated short wave radiation on a horizontal surface at the 
elevation and latitude of the snow area in langleys per day. It may be 
estimated from observed station values near the area or indirectly 
from duration of sunshine, cloud cover, or diurnal temperature varia- 
tion. 


7 is the estimated albedo of the snow surface. Values of a can be esti- 
mated from figure 1. A melting snow has a smalier albedo than the 
non-melting snow; the albedo of new dry snow may be as high as 0.85 
and falls very rapidly to lower values with time, following a decay- 
type curve. The accumulation of residual impurities, as snow melts, 
tends to decrease the albedo and in forests it. may be as low as 0.35. 


¥ is an estimated average basin cover, effective in shading the area 
from solar radiation, estimated as a decimal fraction. The estimate 
of F involves the consideration of the type, height and density of the 
forest cover. 


T is a factor intended to correct I; for an average slope angle and 
orientation of the snow covered area. The effective short wave radia- 
tion decreases as a sloped surface is oriented from south toward 
north, and as the slope angle is farther away from normal to the sun’s 
rays. In effect k! = I/lj, the ratio of the insolation on the true snow 
surface to that on its horizontal projection. Figure 2 shows I/Ij values 
for the entire year for a southerly and for a northerly oriented unit 
area of 25° slopes in 46° 30 N Latitude at 5200 feet above mean sea 
level. 


is an average corrective factor for the degree of exposure of the snow 
area to wind. This factor has a maximum value at 1.0, but it may be 
as low as 0.3 for a forested area. 
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MEAN DAILY ALBEDO, PERCENT 


8 
AGE OF SNOW SURFACE, DAYS 


VARIATION IN ALBEDO WITH TIME 


FIGURE | 


N is the estimated cloud cover, expressed as a decimal fraction. 


Attention is called to equation 19 in which k', Ij, v, and k are missing. In 
a heavily forested area, the effect of wind is dampened to a near constant 
value and insolation is a small portion of the total heat supply, approximately 
equal to evapotranspiration. The heat energy, effecting the generated snow- 
melt runoff, essentially includes convection-condensation and long wave radi: 
tion. These components are satisfactorily defined by a linear function of the 
air temperature (T) and the dewpoint (Tg), within the range experienced dur- 
ing the melt period. A combination of both the rational and statistical ap- 
proach resulted in the coefficients shown in equation (19) for the generated 
daily snowmelt runoff. 

Equation (20), for the forested area, has been derived by a rational ap- 
proach, assuming insolation melt plus ground melt equal to evapotranspiratic 
and therefore, the generated daily snowmelt runoff (Mg) equals the convectio! 
condensation melt plus the longwave melt. 

When equations (19) and (20) are used, “water excess” for routing is ob- 
tained by subtracting only the infiltration water which is not available to dire 
runoff within the hydrograph period. Equations (21) and (22) are derived by 
rational approach, considering point surface melt. The snowmelt M is the 
water input at the surface and includes the evapotranspiration. 

The snowmelt during rain was derived from clear weather melt equations 
Since Tg = T, is a permissible substitution during a rain-on-snow event, the 


daily cangee! @niseneansation melt from equation (18) for an open to partly 
forested area is: 


1s = 0.0084 k Vs0 (Gs = 32) (2 


ASCE SNOW HYDROLOGY {al 


w 


w ia TENF 
< | | 
ra + Nore. + —— paieilhee 
& a +— : = — 
a Curves give theoretical! volues | [ Direct solar radiation only | 
2 | caleuleted for 46°30'N latitude |_|—-——| Direct solor plus diffuse | 
x ond for on elevation Of approx. | sky radiation. | 
z 5200 F4. MSL (of ter Hoeck). 
N A | 
3 
: | | | | | 
< | I] 
> | ‘6 
ae T aoa 
Zz =~ South s/ope, 245° | 
e oN SF 4° ! t- AL 
< | < | 
a S 

+— = —} 4 
2 — we a 
fo) x | ua 
RNS a en ane 
E i | 
—_ + = j : 
: = ato 
a1 Horizontal, Bost or West s/opes=/ | 
2 | 
aw + | 
a —— = — 
° Ye Ss ~~ a! 
= = ~ Of 
g 2 a 
z a >| 

ot <4 ? 

< : ai 
a 39 N oat ees 
2 al y 
2 | L| 


DAILY SOLAR RADIATION ON CLEAR DAYS ON NORTH AND 
SOUTH SLOPES OF 25-DEGREE GRADIENT RELATIVE TO 
RADIATION ON A HORIZONTAL SURFACE 


FIGURE 2 


For a forested to heavily forested area, where wind may be considered a con- 


stant, a statistical treatment gave 
M' = 0.045 (T. - 32) (24) 
ce a 


Assuming the temperature of the cloud base, forest foliage, and the air 
above the snow are the same, the long wave radiation melt would be 
mre a (25) 
M' | = 0.029 (7 - 32) 
As a result of rain heat, assuming the air temperature = the rain tempera- 
ture, the rain melt y 


M'. = 0.007 (T, ~ 32) P. (26) 


where P,. is the daily rainfall. 
The short wave melt varies with the intensity of the rain. But it is very 


small and it is considered satisfactory to express it as a constant or as a 
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function of the forest cover factor F, like 

M';g = 0.05 for forested to heavily forested areas. 

M!l ag = 0.10 for lightly forested and open areas, or 

M'rg = 0.05 (2-F) (27) 


Ground melt is considered negligible, if any, during a time when the ground 
heat supply may be used up by the moving surface and subsurface water. Melt 
by rain impact (kinetic energy) is insignificant. 

The components are combined, and the following equations are written to 
represent the daily snowmelt values during a rain-on-snow event: 

For open to lightly forested areas, using equations (23), (25), (26) and (27) 


M' = (T, - 32)(0.0084 kVs5q 4 0.007 P. # 0.029) # 0.05 (2-F) (28) 


For forested to heavily forested areas, using equations (24), (25), (26) and 
(27) 


M' = (T, -32)(0.074 # 0.007 P.) # 0.05 (2-F) (29) 


Note that for a given basin k and F are constants. When they are determined, 
the equations are further simplified, allowing the construction of a simple 
table or nomograph for zonal or basin melt values. 
The snowmelt over the entire basin (Mp) is: 
M, = MA, 
A, 


where Ag/ A, iS the ratio of the snow covered area to the area of the entire 
basin. b 


(30) 


Effect of Snow Pack Conditions on Runoff 


Water passing through a snowpack may freeze if the snow is “cold” and 
dry; it may be retained against gravity as hygroscopic and capillary water; 
it may be retarded or detained by ice dams, by choking, and by suspension in 
a form of slush; or it may continue in well-defined channels. In a mountainou 
area of snow cover all of these phenomena may occur simultaneously during 
winter and early spring. 

As in other mediums, water will move through the snowpack in a path of 
least resistance. The rate of movement and the storage of water in snow 
depends on the topography, the temperature, depth, density and the crystaline 
structure of the snowpack. Regardless of its path, before water can pene- 
trate through a layer of the snowpack, the path in that layer is brought to 0°C 
by the heat of fusion of the freezing water; each snow particle in the path is 
coated with water, and the capillary spaces filled. Only the excess continues 
downward. Upon reaching an impervious ice plane, the water flows over it 
until a weakness in the ice plane allows part or all to pass through to the next 
snow layer below. Cells of cold snow between ice planes may still remain in 
the pack when water reaches the ground level. However, continued rain and 
melt will eventually disintegrate all ice planes, change the crystaline struc- 
ture and render the entire pack isothermal and saturated to capacity. At this 
stage, the water will move downward almost vertically, and the snowpack will 
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ield all input to runoff with little delay. Water freezing within the pack and 
ter held against gravity (not including any possible slush water in level 
eas) constitute the retention storage or the initial loss to direct runoff. 
ch water is retrieved only as the snow melts. The magnitude can be com- 
puted, and the time delay determined from the input intensity. 
The cold content is the heat required to raise the temperature of a “cold” 
nowpack to 0°C. Its water equivalent is the amount of water that would sup- 
ply the required heat (80 calories per gram) by freezing within the pack. 
Taking the specific heat of snow at 0.5, the cold content (h) in calories is: 


h = 1.27WT; (31) 
os 
and its water equivalent (W,) in inches is: 


W. = 0.00625 WT, (32) 


when W, is the initial snowpack water equivalent in inches and Tg, is the 
average snow temperature in degrees below 0°C. If a known amount of water 
of 0°C entered a “cold” snowpack of uniform texture and density, the gen- 
erated runoff appearing at the bottom of the pack would be less than the inflow 
by an amount equal to the water equivalent of the cold content (W,) plus the 
hygroscopic and capillary capacity of the snow. 

Like soil, the snowpack is capable of supporting hygroscopic and capillary 
water in amount varying with the mass of the pack, which is a function of 
size, Spacing and shape of the snow crystals. The density, as an integrating 
factor of the structure of the snowpack, appears to be an index of the amount 
of such liquid water. Laboratory experiments with melting snowpack of 0.43 
(0 0.56 g/cc, after complete drainage, showed retention of liquid water 
amounts, held against gravity, from 2 percent to 7 percent by weight. No 
studies are known to show the order of magnitude of the free water holding 
sapacity of the winter snowpack which has generally an average density less 
han 0.43 g/cc. The light winter snowpack with small grains has greater sur- 
lace area to be wetted but larger voids which may not be capable of supporting 
sapillary water. Light dry snow, of course, can not remain light after drain- 
we. When rain falls on it, the snow immediately settles and reaches a den- 
sity near 0.25 g/cc after which the settlement rate is relatively slow. Since 
here is a lack of information on the free water holding capacity of the winter 
snowpack of relatively light density, the lower limit (3 percent) of the experi- 
nental results is considered satisfactory. In general then the liquid water 
lolding capacity of the snowpack against gravity (after drainage) is stated as 


Spf", (Wf W,) (33) 


n which S¢ is the water stored and f!!,, is the liquid-water-holding capacity 
f the snowpack in percent of W. In the winter, a value of fl", = 0.03 should 
rove satisfactory, and in the spring for the daily thawing of the night crust 

| value of f ", = 0.04 or 0.05 is recommended. 

_ The magnitude of the retention storage (Sy), as an initial loss, may now be 
‘epresented as the sum of the water equivalent of the cold content and the 
iquid water holding capacity of the snowpack, or 


S, = 0.00625 WT, 4 £", (Wot W) (34) 
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In a natural ripe snowpack in which the planes are disintegrated and the 
pack completely primed to yield runoff, the water travels vertically at about 
30 to 60 inches per hour depending inversely on the density; the denser the 
snow, the slower the velocity. But the time of travel (tj) in a ripe snow is 


t;-= D/V (3 


in which D is the depth and V is the velocity of water in the snowpack. The 
velocity of water in a ripe snowpack is relatively large compared to the 
velocity in a cold and dry winter snow. In a cold snowpack the time of travel 
(t) depends on the intensity of the inflow or rainfall (r) plus melt (m), the re- 
tention storage capacity, and depth of the pack. That is, 


r #D 
c7#m V (3 


Where Sy, is the amount shown in equation (34), and t is the total time elapse 
from the beginning of rain to the time of appearance of appreciable runoff at 
the bottom of the pack. The relative magnitude of the retention storage, (Sy) 
and time delay (t) are shown by the following examples: 


(1) Given: a cold winter snowpack at 7,000 ft. elevation where the snow 
depth is 50 inches with an average density of 0.30 g/cc and temperature of 
-5°C; rainfall and melt of 2.40 inches is forecast for the 24-hour period at 
an average rate of 0.1 inch per hour. Find Wo, the snowpack water equiva- 
lent, Sy, and t. 


Wo = 50 (0.30) = 15 


Using f'', = 0.03 and V = 40 ins. per hour, S, = 0.00625 (15) 5 + 0.03 / (15) 
0.00625 (15) 5/7 = 0.93 ins., leaving 1.47 inches for runoff which will begin 
0.93 50 


to appear under the pack at this site after t = 0.10 + 40 


= 10.5 hours from tl 
beginning of the rain. 
At time of appearance of water under the pack, approximately 0.13 in. 


(=79 (0.10)) is in the snowpack, moving downward. If this was the first rain 


the snowpack had experienced, t >10.5 because of its detention storage capa 


ty in the form of several inches of slush at the bottom of the pack and perha) 
on several ice planes. 


(2) Given: At an open site on a clear spring night, the top 6 inches ofa , 
snowpack cooled progressively to an average temperature of -3°C by 0500 
hours (sunrise); this cold layer and the remainder of the pack had a density 
of 0.50 g/cc; below this crust the pack was drained of water except for its 
retention storage. From solar radiation and temperature forecast, it is est 
mated that the melting will begin at 0530 hours, its intensity varying with ti 
between 0530 and 1100 hours, as follows: m =0.01 t2 in/hr. Find S,. Usin 
p = 0.05, the night crust will require Sy = 0.00625 (6) (0.5) + 0.04 (6) (0.5) = 
0.18 in. It can be seen that ripening the snow crust each day requires an ap 
preciable amount of the daily melt. 

Note that the average temperature and the water equivalent of the snow- 


pack must be known or estimated, before the retention storage and time del: 
can be computed. 
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Application of Principles within Columbia Basin 


The principles and procedures developed as a result of the Snow Investiga- 
ions have assisted in solving problems in snow hydrology in those areas of 

€ mountainous western states investigated by the Corps of Engineers. Sud- 
en melting of a moderate snow cover over the Northern Plains areas can 
ause excessive floods, and some of the actions and reactions in these cases 

e not yet completely solved. In the Columbia Basin, use of these principles 

assisted in obtaining seasonal volume runoff forecast procedures, stand- 
id project floods, spillway design floods, and snowmelt relationships and in- 
lexes. 

The Detroit project, located on a tributary to the Willamette River, is a 
nulti-purpose project in which runoff from snowmelt floods is detained and 
ised for the generation of power. The North Santiam above Detroit Reservoir 
s located on the west slope of the Cascade Mountains about 60 miles south- 
ast of Portland, Oregon, and has a drainage area of 438 square miles. Ele- 
ations range from 1200 to 10,495 feet above mean sea level, the mean eleva- 
ion being 3718 feet. A large percentage of the area is comprised of valleys 
nd ridges with steep slopes, and a heavy stand of coniferous trees covers 
nost of the area. Soil cover is relatively thin, but there is considerable duff 
nd litter under the heavy forest canopy. The climate is characterized by wet, 
noderately cold winters and dry warm summers. Normal annual precipita- 
ion over the area is estimated to be 82 inches, ranging from below 70 near 
he Detroit project to well over 100 at the crest, where snow accumulates to 
reat depths. About 90 percent of the precipitation occurs during the October 
hrough May period. 

The water balance method, discussed earlier in this paper, was used for 
orecasting the volume of runoff from the basin above this project, and the 
rocedure used is considered representative of a method adaptable to an 
verage project basin. Several of the steps taken in this procedure deviate 
rom the procedure derived for a laboratory area. No estimate was made 
Ir losses by interception, since these are smail compared to the total rain- 
ull, which in turn is not accurately known, having been estimated by summing 
1e generated runoff and the calculated evapotranspiration loss. Likewise, no 
stimate was made of gage-catch deficiency due-to the wind. 

Precipitation, snowfall, streamflow, air temperature, and snowpack water 
quivalent data are available for varying periods. The only adequate temper- 
ture record available is that at Detroit, necessitating the use of lapse rates 
) obtain estimated temperatures at higher levels. Precipitation data are 
vailable for five stations, one of which was begun in 1909. Water equivalent 
; measured at four snow courses having records since 1941. Because of 
culation during the construction of Detroit Dam, the streamflow record 
ibsequent to 1951 was not considered usable for this investigation, and the 
sriod of record suitable for study extends from the beginning of snow-course 
scords (1941) through 1951. 

The snowpack remaining on the basin as of August 31 each year is negligi- 
e, and the reservoir is filled after the threat of a rain-type winter flood has 
rgely passed, so a reliable forecast procedure for the runoff available from 
ebruary 1 through August 31 was needed. Forecasts covering the periods 
arch-August and April-August were also needed. Equation (4) is used for 

l periods. The basic equation used for all periods is as follows: 
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Q:=P# (W, - Wo) #(Q, - @) - L (3° 
in which Q is the available runoff during the period from the date of the fore 
cast to August 31. P is precipitation, Wi and W2 are the initial and final 
snowpack water equivalents respectively, Q; and Q9 are the runoff volumes 
under the initial and terminal recession curves, and L is loss, all values 
being estimated as inches of depth over the drainage area. The final snow- 
pack water equivalent, Wy is equal to zero in this case. The snowpack 


water equivalent for the date of the forecast (W1) is obtained from a plot of 
the water-equivalent at each of three key stations, at known elevations, from 
which an indication of the water equivalent for each 10 percent of the basin 
area is read. The average snowpack water equivalent is obtained by appli- 
cation of a factor derived from a study of the relationship of precipitation, 
loss and runoff during the 11-year period of study. Q, is the recession 
volume corresponding to the discharge on the date of the forecast and Qg is 
the recession volume corresponding to the discharge on the last day of the 
forecast period. Losses, L, were computed by Thornthwaite’s method for 
each year, based on temperatures at Detroit. This method relates losses to 
climate, and estimates losses by an empirical formula using only air tem- 
perature as a variable, with adjustments for precipitation experienced and 
for the latitude of the area considered, which indicates the length of the day. 

Precipitation during the forecast period must be accounted for in use of 
the forecasting procedure, since spring precipitation constitutes a major po! 
tion of the water supply during this period. Values of precipitation may be 
assumed on the basis of normal or critical conditions during the spring peri 
or someday it may be possible to utilize long-range precipitation forecasts. 

On February 1, the initial snowpack water equivalent W, amounts to only 
about 25 percent of the total of Wy plus P; this forecast is greatly dependent 
upon future climatological events and is generally to be considered a prelim 
nary estimate. On April 1, the value of Wy is about 60 percent of the total o 
Wj plus P, and the forecast of this date is reliable enough to be of major as: 
sistance in the operation and refilling of the Detroit reservoir. 

The following is the February 1-August 31, 1958 runoff forecast of the 
water balance method: 


February 1, 1958 (11-year average - 13.2 ins.), Wy, = + 14.6 in. 


February 1-August 31 normal precipitation, P =+ 30.6 in. 
February 1-August 31 normal losses, L =- 8.8 in. 
Recession runoff volume for observed February 1 

discharge (6227 cfs) Q; =+ 9.2 in. 
Recession runoff volume for normal August 31 

discharge (633 cfs) Q2g =- 5.1 in. 
February 1-August 31 runoff forecast, Q = 40.5 in. 


Cougar Dam, now under construction on a tributary to the Willamette 
River, is to be a multi-purpose project storing water for flood control and 
power production. A standard project flood is derived for projects similar 
to this, in order to determine the amount of flood control storage that shouk 
be provided if economically feasible. This flood is one which would be 
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exceeded only upon rare occasions caused by a rainstorm falling upon a snow 
covered basin. It requires the combination of rather severe antecedent con- 
ditions and a temperature sequence that would cause some of the snow to melt 
and contribute to the flood. The largest flood recorded at the dam site oc- 
curred in 1953, and this rain-on-snow event was considered to precede the 
standard project flood. 

The 210 square-mile drainage area above the Cougar dam site, elevation 
1532 feet mean sea level was divided into five elevation bands of about 1000 
|feet each. Analysis of snow surveys and the climate of the region indicated 
‘that the basin could reasonably be expected to have a 100 percent snow cover 
with 8.2 inches of water content at the dam site and an increase of 7.5 inches 
of water content for each of the five bands. This is considerably more than 
the normal and could accumulate only under temperatures less than normal. 
It was conservatively assumed that the basin snow cover was at 32°F and re- 
quired only hygroscopic and capillary water of about 2 percent by weight be- 
fore yielding runoff, amounting to an initial loss of 0.55 inches, part of which 
was regained as the snow melted. Melt due to short wave radiation and 
ground heat is insignificant during a rainstorm and was omitted. The snow- 
melt was computed by zones, using a melt rate of 0.08 inches per degree-day 
above 32°F applied to appropriate air temperatures for each zone. The melt 
rate conforms to that previously determined at the Willamette Basin Snow 
Laboratory for the condition of rain-on-snow in heavily forested areas. Snow- 
melt from rain itself was added separately. The hydrograph of this flood 
would have a peak discharge of 40,000 c.f.s. or 190 c.f.s. per square mile, 
and a total volume of 194,000 acre-feet. Storage of 154,000 acre-feet of run- 
off would be required to control this flood, and a usable capacity of 155,000 
acre-feet is being provided. 

Some of the variables for clear-sky snowmelt are neglected in this compu- 
tation. The heavy forest on this basin markedly reduces the wind speed varia- 
tion over the snow surface. Therefore, a wind variable is unnecessary, and 
evaluation of convective melt can be accomplished by use of air temperatures 
alone. With intense rainfall the air would be saturated, so that the dewpoint 
temperature would be equal to the air temperature. For this reason, conden- 
sation melt may also be evaluated from an air temperature function alone. 
Rainfall over the basin during the most intense 6-hour period was estimated 
to be 3.26 inches, and snow melt during this same period was 0.72 inches. 
Total rainfall for the design storm amounted to 10.75 inches, and snowmelt 
totalled 3.64 inches. The temperature at the 3900 foot level was assumed to 
reach a maximum of 49°F. 

Several spillway design floods have been computed for projects in the Col- 
umbia River Basin using the procedures and relationships that have been de- 
veloped as a result of the Snow Investigations. Such a flood was developed for 
the Libby project, on the Kootenai River in Canada, and just recently for 
projects on the Salmon River in Idaho. These determinations require a con- 
siderable amount of detailed work, and they will not be described here. 

A very important application of snow hydrology in the Columbia Basin is 
evolving in the day-by-day use of data gathered by electronic means which 
must be interpreted and used in producing a forecast of flood flows and gage 
heights in the ensuing short-term period. Reservoir regulation for flood 
control and flood fighting operations at local critical areas and in the reach 
of leveed river downstream from Bonneville Dam will be assisted thereby. 
one a forecast must be completed within a short time after the raw data 
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become available to be of use, the hand or desk machine computations enter- 
ing into them have been restricted to a relatively few. With the advent of a 
high speed digital computer, and a routine written for the Columbia Basin, 
forecasts will be made for a large number of small tributary basins, based 
upon antecedent weather and forecasts of temperature and precipitation to be 
expected in the next several days. The same routine will route the flows 
from the small tributaries through lakes and channel storage and will furnish 
plotted hydrographs at key points. Such a routine requires the derivation of 
indexes and coefficients pertaining to each of the small tributary basins. 

One of the tributary basins selected for analysis was that of the Boise 
River above the Twin Springs gaging station, where the 1954 and 1955 spring 
snowmelt hydrographs were used. The Boise River above Twin Springs, 
Idaho drains an area of 830 square miles of the Sawtooth Mountains. The 
area has deep valleys, steep slopes and narrow sharp-topped ridges, ranging 
in elevation from about 3000 to 9000 feet above mean sea level. The stream 
system and instrumentation are shown in Fig. 3. Forest cover consists 
principally of conifers which are estimated to have a shading effect on 30 per 
cent of the basin area. Snowmelts for each day of the periods were computed 
by the thermal budget index and the flows resulting from this released water 
plus the precipitation for that day were routed to the gaging station for com- 
parison with the observed hydrograph. This comparison for the May-June 
1955 period is also shown in Fig. 3. The index derived for the snowmelt ove1 
this basin is: 


M = 0.0267 T| 4 0.002276 - 2.00 (38 


where M is the daily snowmelt in inches over the basin, Tmax is the daily 
maximum temperature, at Idaho City, in degrees F, and G is a radiation 
parameter, in langleys (equal to absorbed shortwave radiation plus esti- 
mated longwave loss.) 

It will be noted that one of the terms in the above equation requires that 
a forecast of the radiation to be experienced on ensuing days must be availa- 
ble. Radiation is measured daily at Boise, Idaho and these data were used in 
deriving the equation, but large parts of the Columbia Basin do not have such 
measurements available. This refinement may be used sometime in the fu- 
ture, when more reliable forecasts of the temperature and precipitation and 
other weather factors are possible and when data on radiation become more 
widespread. For the present, a simple maximum daily temperature index of 
snowmelt is being used, which provides fairly reliable estimates of snowmelt 
runoff for day-to-day forecasts of streamflow. 

Another application of indexes to the predication of flood volumes resultec 
from the Snow Investigations. The method is based primarily upon the rela- 
tionship between winter runoff of low-elevation drainage basins in western 
Washington and Oregon, and the spring snowmelt runoff of the Columbia 
River at The Dalles. Indexes of winter temperature and spring precipitation 
are included in the forecast procedures as secondary parameters. The use 
of low-elevation winter flow as an index is confined to regions where the low 
elevation and the high-elevation areas have a common source of moisture. 
Such a situation exists in the region composed of the Columbia Basin and 
western Washington and Oregon, the entire region being well centered in the 
belt of prevailing westerlies. Moisture is carried in a generally eastward | 
direction from the Pacific Ocean, and the amount deposited is largely 
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dependent upon the rate of flow and the preciptable moisture content of the 
air, reflected by both winter streamflow from low elevation areas and the 
accumulation of snow at high areas. Several tests indicated that a correla- 
tion of winter flow from two streams each in Washington and Oregon would 
yield results as accurate with regard to historical data as those which use 
precipitation and snow course data for the principal indexes. This method 
does not produce any information as to the distribution of the snow pack over 
the smaller interior basins of the Columbia River above The Dalles, Oregon, 
but it does produce an acceptable forecast for the basin as a whole quite earl} 
in the accumulation period. 


SUMMARY AND CONCLUSIONS 


The significant developments in applied snow hydrology which have re- 
sulted from the work of Snow Investigations research and analysis are: 


A point snowmelt equation based on heat supply from short and longwave 
radiation, convection and condensation. 


A general snowmelt formula applicable to any drainage area when overall 
coefficients characteristic of that basin are evaluated for the degree of ex- 
posure to wind and radiation. 


Relative magnitudes of melt components under different sky and forest 
cover conditions. 


Numerical appraisal of the storage and delay effect of the snow pack on 
runoff. 


Runoff forecast procedures based on water-balance and on a coastal 
winter-flow index. 


In addition to the above direct applications, numerous investigations were 
performed which were of a more basic nature. These deal with the theoreti- 
cal aspects of snow physics, precipitation, snow cover, and soil moisture 
measurements, forest effects on snow accumulation and evapotranspiration 
loss, heat exchange between the snowpack and its environment, the distribu- 
tion of the snowpack in mountainous areas, and the water balance in areas of 
snow accumulation. 

Our knowledge of snowmelt is not complete, but we know enough now to 
supplement the degree-day approach by including humidity, wind speed, and 
radiation values even if they must be partly estimated, particularly for clear 
weather melt in open and partly forested areas. Only during a rain ina 
dense forest, where the wind speed and the solar radiation are small and 
nearly constant, can the temperature alone define the snowmelt properly. 

Direct application of these principles is now being made in estimating de- 
sign floods for determining the capacity of flood control storage projects and 
river channels, and the related problems of conduit capacity, revetment de-_ 
sign, levee heights, and day-to-day operation of the reservoirs. Seasonal 
operation is based on results of monthly forecasts of the remaining runoff by 
the water balance method. 

The network of stations is being expanded to include more hydrologic in- 
formation and better areal coverage. Radiation data are observed at only a 
few stations and this class of observations should be extended to cover more 
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area. The snow fields in the high mountains are uninhabited, and observa- 
tions of climatoligical events in this large and important area are difficult to 
obtain but steady advances are being made by use of radio-reporting facili- 
ties. The distribution of snow cover during the melt period is also being ob- 
served on a systematic basis by means of aerial snow reconnaissance. 
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FLOW THROUGH CIRCULAR WEIRS® 
Closure by J. C. Stevens 


J. C. STEVENS,! M. ASCE.—The author is gratified that six engineers 
fell versed in mathematics and hydraulics took time to send in discussions 
f this paper. 

The experiments by Mr. Manabar are very interesting even though his 
ipe was small. The sharp edge of the inlet of a pipe on a steep slope is not 
1€ Same as a Sharp-edged circular weir through which free flow obtains. 
/hen partly full the inflow was partially submerged and this is correctly as- 
igned by him as the reason his average flow coefficient of .52 was less than 
1e author’s average of .59. His Fig. 2 indicates that the stage-discharge re- 
ution is virtually a straight line for higher values of H/D. This is in con- 
rmity with the graphs for the laboratory flows listed in Appendix A. 

The head on circular weirs need not be limited at all to H<D for the weir 
lay be surcharged; also it may be submerged and it is often necessary to 
alculate the flow under these conditions. 

For surcharged circular weirs there are two conditions to be met. In the 
ange H<D flow is readily calculated from the elliptic formula with the aid 
[ the tables in the author’s paper. When H>D the weir becomes an orifice 
ir which the common orifice formula is applicable, the head being measured 
‘om the center. Fig. 3 shows the calibration of a §-inch sharp-edged circu- 
1 weir by E. R. Dodge to a head of H=D (Appendix A). However, his thesis(3) 
ives the flow for surcharged heads. The weir was in the outlet end of and 
mcentric with a 15-inch pipe in which the approach velocity was negligible. 
he average flow coefficient of .607 given in Appendix A for this weir was 
sed with the Equation (3) to define the solid line Fig. 3 (a) for H>D. For 
ircharged heads the same coefficient was used in the orifice formula viz 
= .607 AJ2gh where A = .197 sq ft and h is the head on the center of the 
rcular weir-orifice. The small circles represent Dodge’s laboratory flows. 

‘Mr. McPherson has given the effect of surcharging the 1-foot weir from 
mori’s experiments in his Table B. Fig. 4 is a graph of the flow through 

e Omori weir calculated by Eq (3) up to H/D = 1.00 using the average flow 
efficients of the last column in Table A. For the surcharged head the cor- 
‘sponding average coefficient for the weir, .596, was also used. Note the 
noothness of the transition of flow weir to orifice in both examples which 
e well substantiated by laboratory experiments. 

For submerged flow the Villemonte formula (4) may be used. For a weir 
any kind subject to submersion in a stable channel observations will soon 
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define a relation between headwater and tailwater. To illustrate, Fig. 3 (b) 
shows such a relationship (similar to one found for the outlet of a grit chan- 
nel with a proportional sharp edged weir at its outlet.) 

The general form of the Villemonte Formula is 


m 
-Q_ = k(1 - Qb) ( 
Q, vote 


Where Q is the observed flow -cfs. 
Qa, the flow corresponding to the headwater (above) 
Qp, the flow corresponding to the tailwater (below) 
k, a constant 
m, an exponent to be determined 


Both Q, and Q) are taken from the free flow curve of Fig. 3 (a). For most 
sharp-edged weirs k is practically 1.0 and m has the average value of .385 
found by Dr. Villemonte’s laboratory experiments. A nomogram is readily 
constructed for the Villemonte formula for any set of stable conditions.(15) 

In Fig. 3 the dashed line shows the effect of submerging the weir-orifice. 
This well illustrates the versatility of the circular weir-orifice. What other 
type of open channel measuring device lends itself to reasonable accuracy 
throughout a head range from zero up to a hundred or more feet? Even the 
Venturi tube is not accurate for low rates of flow. 

The author is grateful indeed for the errata in Mr. Blaisdell’s discussion 
As one grows older his facility for errors increases markedly, especially 
with the slide rule. The flow experiments have been checked and revised 
tables of Appendix A, where errors were found, are appended hereto. 

The table of Appendix C is revised in its entirety in order to eliminate th 
minor errors in the 4th decimal place for which Mr. Moody takes the author 
to task and quite justifiably so. Any table should be correct to the last deci- 
mal. Mr. Moody’s admirable table of values of Ato 4 decimals with consis- 
tent 2nd is a step in the right direction but interpolation may be necessary. 
The author believes that k? should be carried to 3 decimal places. 

Mr. Kolupaila’s table gives the theoretical flow directly for a 1-foot diar 
ter weir. The argument x has been substituted for k2 in the author’s Eq (3 
The theoretical flow through a circular weir of any other diameter D is ob- 
tained by multiplying the tabular value by p5/2, 

In answer to Mr. Blaisdell and doubtless others who noticed the lack of i 
formation regarding the approach conditions for the weir data in Appendix A 
references for the experiments now include this information. 

The author knew nothing about the Byrd and Freidman Handbook of Ellipt 
Integrals for Engineers and Physicists, Moody’s Ref (1). He has since pur- 
chased a copy and was quite non-plussed to find included 2 errata pages of 
fine print — about a hundred corrections. In the sweet-bye-and-bye perhap: 
there’ll be no errors - “Tis a consumation devoutly to be wished.” 

In weir computation the orthodox method followed throughout the entire 
history of hydraulics has been to compute the theoretical flows, then compa 
them with observed discharges and correct the computed flows by a factor 
which might be called a coefficient of ignorance. 

For all sharp edged weirs that coefficient is some measure of the effect 
the sharp edges distorting the stream lines as water passes the weir. Mr. 
Blaisdell’s use of his formula (a) for circular weirs is unorthodox but cer- 
tainly permissible. To be utterly capricious one might even measure H fro 
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y arbitrary datum, reduce or increase a by a constant or a percentage, 

en use the value of g from Mars or Venus and obtain “goofy” discharges 
nuts) Computed from (A). Parallel observed discharges would then deter - 
ine a series of coefficients that would define a “goofy” curve but from which 
e true discharges Q could doubtless be found. 

The average coefficient C = .365 (Fig. A), he admits gets top-heavy beyond 
/D = .8. Dividing the scattered values of his c by the author’s correspond- 
g scattered values of C (Fig. 2) he obtains the tantalizing curve of Fig. B 
hich the author at first mistook for an improved set of discharge coefficients 
plicable to circular weirs. However, Mr. Blaisdell, by letter, disclaims 
my such intention and states “the values in my Fig. B are a measure of the 
acorrectness of my approach, part of which I suggested might result from 
ay neglect of what Professor Mavis calls ‘shape factor.’ I worked it out and 
ompared my results with yours as a matter of interest. . . I thought others 
aight also be interested.” (See his Ref. 3). 

The data on discharge coefficients presented by Messrs. McPherson and 
olupaila are shown in Fig. 5. Curve I indicates the average coefficients of 
he 4 experimenters given in Table A of Mr. McPherson’s discussion viz 
rmori, Staus, Staus & Sander, and Greve. 

Curve II gives Dr. Kolupaila’s coefficient formula which he slightly modi- 
ied from that of Staus as follows: 


Cy = 0.570 4 1 + 0.041x 
100x 


here x = H/D. The Staus formula used 0.555 in place of 0.570. 

The author’s constant coefficient of .59 fits Curve I fairly well for moder- 
te to high values of H/D. The deviations are -9% at H/D = .10, -5% at .20, 
nd practically zero above .30. 

Mr. Moody has offered a new and valuable approach to orifice flow by his 
rmulas (3) and (4). In order to avoid confusion it would be preferable to use 
»/H)2 rather than k4. For theoretical flow through vertical circular orifices 


e proposes 
Oa? areas AV2gn (4) 


1577 (D/H)* 
hich readily reduces to 


QO. 21a Ane? (21) 


here Ais the elliptic function given in Appendix C using D/H for the argu- 
ent instead of k2 = H/D, while H is measured from the invert of the orifice 
3 in the weir instead of from its center. Obviously D/H < 1.00. 

Fig. 3 shows by a dotted line the results of using (21) for the flow through 
odge’s 6-inch orifice, using the average weir coefficient of 607. (Appendix 
). Fig. 4 shows the corresponding flow for Omori’s 1-foot orifice in Table B 
sing the average coefficient for Omori’s weir experiments of .583 given in 
e first column of Table A. It appears therefore that the average discharge 
efficient found for the circular weir may be used as the average coefficient 
vr the same weir when it is surcharged becoming an orifice as is indicated 
7 the four examples herein given. More laboratory aperiments on sur- 
larged circular weirs would be fruitful. 
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Pages 1455-5&6 - Substitute the following for the corresponding numbered 
references: 


il 


Cone, V. M. Circular Weirs. Journal of Agricultural Research USDA 
March 6, 1916. The approach flume was 20 feet long 10 feet wide and 6 
feet deep. For the 2-foot weir the flume bottom was 4-1/2 feet below the 
center of the weir. 


. Greve, F. V., Flow of Water through Circular Parabolic and Triangular 


Vertical Notch Weirs, Engineering Bulletin, Research Series No. 40. 
Purdue University, Lafayette, Indiana. This bulletin also contains a 
bibliography with summaries but no mention of European experiments. 
The approach channel was 16 feet long and 5 feet wide. The depth below 
the weir is not given but it scales 2-1/2 feet below the invert of the large 
weir. 


. Dodge, Eldon R. A Thesis Submitted for a Masters Degree, University of 


Wisconsin, Madison, Wisconsin. 1935. (Not published). The sharp- 
edged circular weirs reported herein were concentrically in the end of tw 
sizes of tile sewer pipe as follows: 


Weir Nominal D Pipe I D Weir Nominal D Pipe ID 
2-1/2 -inch 8-inch 4-inch 15-inch 


4 -inch 8-inch 6-inch 15-inch 


. Thijsse, J. Th. Discharge Meters Director Hydraulic Laboratory 


(Waterloopkundig Laboratorium) Delft, The Netherlands. Three sizes of 
flumes were used. Weirs were not rated beyond H/D = .75. For the 
2.36-foot (72 cm) weir the nominal dimensions of the flume were, width 3 
feet, depth below invert of weir 1.0 ft length of 9 ft. There is strong evi- 
dence that the increase in the coefficient for all 3 sizes is due to velocity 
of approach. However these weirs were calibrated in place. 


For the 1.77 foot (54 cm) weir the flume dimensions were width 2.2 feet, 


depth below invert of the weir .75 feet, length 7.0 feet. 


For the 1.181 foot (36 cm) weir the flume length was 4.6 feet, width 1.5 


feet, depth below invert 0.5 feet. 


Page 6 add the following references: 


14. Submerged-weir Discharge Studies by James R. Villemonte, Engineerin; 


News-Record Dec. 25, 1947. 


15. Proportional Weirs for Sedimentation Tanks, by J. C. Stevens, Proc. 


Paper 1015, June 1956. 


Appendix A - Substitute the following tables for those of the corresponding 
Author and diameter. 
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E. R. Dodge 
H H/D Q Q 
.084 .398 .018 012 
100 .474 -0£5 016 
108 512 .029 019 
a &y-4 .626 042 .026 
. 143 .677 B48 .03.3 
150 far a WE 052 = .031 
-152 -720 .053 032 
.169 .802 -063 .036 
171 810 064 .039 
177 839 -068 .043 
-189 .896 075 .046 
aL9S 915 -077 .050 
.202 -958 .083 .051 
211 1.000 2084 ) 055 

E. R. Dodge 

.128 .380 054 .034 
-142 -423 .066 .040 
-161 478 -082 .050 
.168 -500 4089 . 055 
219 .652 -144 .084 
.229 .682 -155 .097 
235 .700 .162 .096 
.243 es) .174 .101 
.262 -780 »L95- 5194 
RE) .812 5208) .12zZ 
.294 .875 233° .136 
.299 .888 PZage 39) 
.302 -900 244 .151 
pond -955 265° .156 
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.628 
-606 
.610 
.618 
-583 
.625 
-593 
-580 
.584 
-587 
585 
-583 
.619 
-589 
599 


-50 


-50 


-20 
25 
.30 
+35 

-40 
45 
-50 


-40 
-50 
760 
.70 
.80 
-90 
1.00 


F. W. Greve 


.628 
.540 
568 
UA 
.550 
.580 
580 
574 
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F. W. Greve F. W. Greve 
D H H/D Qt Q a D H H/D Qty Q c 
12333 280 E210 9540" 2522. 7.597 si 148 159 2 O98 204 .128 .627 
922, eaa2 .709 417 588 207 .118 345 .214 .620 
‘415° 131) . bok 49 675 588 254 .145 520 .316 .608 
495 .372 1.611 .935 .580 -321 1184 = 825.4820 597 
577 .433 2.139 1.244 .580 388 .222 1.183 .910 .600 
582 .437 2.176 1.260 .580 435.249 1.476 .875 .593 
703 .527 3.069 1.797 .586 487 .279 1.837 1.080 .589 
(814 .610 3.986 2.296  .575 528 .302 2.141 1.250 .584 
.967 .726 5.380 3.070 570 588 .345 2.756 1.540 .558 
1.124 .844 6.873 3.900 .567 -647 .370 3.143 1.830 .582 
1,146. -.860 7077-51995 565 -703 .402 3.6702,140 .583 
1,332 1.000 8.780 5.100 .580 .750 .429 4,142 2.390 .577 
580 801 .458 4.6772.710 .579 
.853 .488 5.257 3.040 .578 
907 .519 5.883 3.400 .579 
590 

V. M. Cone 

LSO0u205- 133 2088 Ae18 .604 2.495.311 .125 .942 .573 .608 
eS) AGT. 46m -—25 535 .409 .164 1.607 .966 .601 
300.200." .Go0mmess 575 453.181 1.943 1.175 .604 
(35 6 7233). B86», 51 575 490 .196 2.265 1.359 .600 
(40 W267 © 1152e 5 {69 598 526 .211 2.608 1.550 .594 
50 .333 1.758 1.03 586 563 .226 2.978 1.765 .593 
60 .400 2.480 1.45 584 585 .235 3.209 1.895 .590 
.70 .467 3,308 1.90 575 620 .248 3.5602.118 .594 
80 .533 4.211 2.44 580 .665 .267 4.105 2.424 .590 
579 .699 .280 4.498 2.648 .589 
723 .290 4.814 2.823 .588 
596 

F. W. Greve V. M. Cone 


2.25 .280 .124 .715 .450 .630 3.00 .30 .100 .953 .48 .504 


.351 .156 1.128 .690 .611 A0-+ 2133) 12686. 928 let 481 
413 .183 1.534 (840 .548 50 .167 2.6401.22 .463 
.525 .233 2.441 1.50 .614 -60 .200 3.7301.72 .462 
.627 .279 3.451 2.02 585 -70 =.233 5.006 2.30 .459 

. 618 80 .267 6.5062.95 .454 


-90 .300 8.142 3.67 .452 
1.00 .333 9.932 4.44 .447 
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F. W. Greve : EF. W. Greve 

D H H/D Qt Q GC Dee ahi) Q Q e 

. . Sao .191 .254 PLO 5 20Y 579 1.00 .205 .205 251 .146 (580 
( =. 23 297 »250) "348 592 254 (254 379, .2202%580 — 

= .308 411 -461 .263 570 2301 630 546 .314 576 

4 ; wale -416 471 .268 .570 364 .364 .754 .430 2577 

: 56. 475. .6032Ne45 572 -430 .430 1.029 .594 .578 | 
.363 484 .624 .356 570 520 .520 1.461 .836 .572 


-575 -526 .526 1.491 .850 .570 
; +634 .634 2.079 1.190 .573 

e132 832 26591520) 572 

-873 .873 3.5302.000 .566 


1.000 1.000 4.2802.520 .590 - ~ 
227 © 194 112.583 576 
269" 269) S158. 35.58%. 3 ; — ae e 
o32 402" 2231 575 Ee 
-348 .440 .254 578 : 
PALS 6205358 —— 5579) aa iVEEM “Conet 7-2 Gitariog 
502  .869 .495 .570 ee. 


590 1.162 .659 575 RS te dele 
i595. 1,183 -.671 ; 
685 1.508 .864 573 
-785 1.894.1.076 569 
794 2928 D2 hae 1S 
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Page 6, 3rd line before Acknowledgments substitute: « 


¥ 
= 
Jorissen, Andre, Le Déversoir Circulaireen Mince Poroi (A Study of the 


Appendix A p 11 Heading of Table for D = 1.00 substitute: 


_ 


Vv. M. Cone 


Appendix A p 12 Heading of Table for D = 1.50 substitute: 


V. M. Cone 


Appendix B p 14 line 2 substitute: 


1 @ = 162 D9/2 (a(n - Rapa) & - 42 - eA) (1 ~ ke) x] ~~ ne 

Appendix C p 17 line 1 substitute a : : ; 

a. Values of A= 2(1 - 2 4 kt) E - (2 - w*) (i ak} Ke ; " , 
A, Fs * ws ™ we 

_ Appendix D p 22 Equation (3) retype 2, : ames So - ea ie 

" p ; , ~ £0 Pa ee ak 

Vie ey sn? x = Dk? sn® x ; ; "ol oe tee Site aie: tank ¢ 

f (os cs _ - Lidys . Sire 4 - 


— ~~ 
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Appendix D p 22 line 4 


y = Hsn* x = Dkésn* x (3) 


Same p22 12th line from bottom substitute: 


1 
= (4p°x8sn4 x cnt x dn* xf 2 gx (9) 


Same p 22 4th line from bottom substitute: 


2 r4 % 
ens = | ~ sn” ,and dna x = l] — sist y 


Same p 24 Eq (20) substitute 


Garam we OY A 
t 75 V2g p [z gt Ey) ee Cee K] (20) 


In Text: 
Page 2 line 16. substitute: 
Experiments by V. M. Cone(1) 


Line 4 from bottom substitute: 
weir flow was V. M. Cone who made his experiments in 1916 at Fort Collins. 


Page 4 line 18 substitute: 
Appendix (C). The following table gives A to 4 decimal places for values of 


k2 to 3. 


Page 5 line 4 substitute: 
The experiments by V. M. Cone are all near the 59 per cent line except 
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0 


0.0000 
0.0003 
0.0012 
0.0026 
0.0047 


0.0073 
0.0104 
0.0142 
0.0185 
0.0233 


0.0287 
0.0346 
0.0411 
0.0481 
0.0557 


0.0637 
0.0723 
0.0814 
0.0910 
0.1011 


0.1117 
0.1228 
0.1344 
0.1465 
0.1590 


0.1721 
0.1856 
0.1995 
0.2140 
0.2288 


0.0000 
0.0004 
0.0013 
0.0028 
0.0049 


0.0076 
0.0108 
0.0146 
0.0189 
0.0238 


0.0293 
0.0353 
0.0418 
0.0489 
0.0564 


0.0646 
0.0732 
0.0823 
0.0920 
0.1022 


0.1128 
0.1240 
0.1356 
0.1477 
0.1603 


0.1734 
0.1870 
0.2010 
0.2154 
0.2303 
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MECHANICAL ANALOGS AID GRAPHICAL FLOOD ROUTING? 


Discussion by C. O. Clark 


nO: CLARK, 4 A.M. ASCE.—The mechanical analogs of aid in graphical 
flood routing presented are of interest to all who must manage the natural 
waters of rivers and streams, either by scheduling and alloting the time of 
availability of beneficial use of that water which is harnessed and under con- 
trol, or in warning of the unavoidable rampages of the public enemy when it 
is not under control. 

The philosophy by which the author searches for his mechanical analogs, 
and finds the dampening effect of storage to be kin to the effect on the move- 
ment of a rear auto wheel of the distance between the wheels, makes refresh- 
ing reading. Such refreshment is needed by those who find the need for think- 
ing and the opportunity for it completely eliminated by the codes, rules and 
regulations, and “accepted” formulas which make up much of the necessary 
pieces of routine engineering work. There are excellent examples of applied 
logic and scientific observation in the manner by which each of these analogs 
has been developed and built into a solid tool. There is food for thought in the 
fact that the author who conceived and built electrical analog computers for 
the solution of these problems, has now turned to mechanical analogs. 

One of the analogs is as small and as subtle as a sliderule, and is as per- 
sonally possessable as a straight edge/or a pencil. It is a truly Euclidian 
tool. Even more, by an as yet unproved modification of this tool, the author 
has opened a promising door toward a variable, non-constant, kind of solution 
of a problem which exists in the nature of streams as they are, but which 
solution has defied the limitations of brains, two-dimensional paper and the 
timeclock up to now, namely the class of problems involving variable relations 
between storage and discharge, variable lag and storage. 

The writer is interested in the template type analogs, but this does not de- 
tract from the more mechanical ones, which are certainly ingenious, probab- 
ly accurate, and based on sound analogy. This interest is very strongly in- 
fluenced by the fact that the underlying principle of the templates, which Mr. 
Kohler now applies to long reaches of the main stem of a river system, is the 
same principle which underlies an application to an entire watershed to de- 
velop the type of flood hydrographs which the writer could call unit hydro- 
graphs, if they had a unit duration of rainfall, but must call instantaneous 
hydrographs when that unit is zero. This essential principle is that two steps 
of simple analysis are first found(1) to be essentially comparable to the more 
complex and simultaneous action of aquatic storage as it occurs in nature and 
as McCarthy approximated it in flood routing formulas developed for the 


Muskingum watersheds in Ohio.(2) 

a. Proc. Paper 1585, April, 1958, by Max A. Kohler. 
1. Hydr. Engr., Tulsa, Okla. 
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The author calls it lag routing, or lagging and then routing through storage, 
which he explains is bodily translation in time of a hydrograph without distor - 
tion, followed by a distortion or attenuation much in the manner effected by 
storage in a large level reservoir where storage is strictly a function of in- 
stantaneous outflow. This is a mental analog, of course, because the one kind 
of a waterway in which none of these mechanical analogs is needed is the level 
reservoir of great depth. Nevertheless, it is a correct analogy and it is the 
limiting case of a class of problems which vary from the headwater problems 
of sewer inlets where translatory effects have been considered large and pond. 
ing effects negligible, down through the streams which embody both in every 
proportion, onward to the large reservoirs, great lakes, and then to the deep 
seas in which the largest waves made by floods are but ripples on a surface 
translated in negligible time and at velocities compared with that of sound, 
and involving no appreciable lateral movement of mass. This idea is not new. 
Translation or lagging was built into flood calculation and thought in the so- 
called rational method of decades ago, repeated but not originated by Ramser 
in 1927. It was applied to urban problems of storm sewer design before the 
turn of the century, one writer says by Mulvaney in 1854, and to the rural 
problems of drainage at the time that Yarnell and Ramser were young men in 
these fields. It is just as natural for a headwater engineer to think solely in 
terms of translation, or a horizontal displacement in time and space, as it is 
for a downstream reservoir or maritime engineer to think of storage and of 
waves as something that involves only vertical displacement in space. Puls, (4) 
working on large tributaries of the Ohio River was one of the first to show 
that even the translation of river waves involved great storage, like reser- 
voirs, but not so effective in attenuation. The writer has tried to knit together 
these upstream ideas of time required for concentration and the downstream 
ideas of storage which must be filled before flow can become steady, into the 
concept of an instantaneous hydrograph(1) which is the least size of hydro- 
graph that could have a peak and the minimum duration if all elements of the 
watershed were equally involved. Hence it is dealing with a minimum volume 
associated with a flow and a duration and an influence of each on the other and 
of upstream and downstream elements on each and both, or the minimum 
wave which is stable enough in form to establish essentially steady flow and 
transmit it along a sloping channel. 

What is definitely new in this work of Kohler’s is the correlation of transla 
tory storage with upstream storage, or lag, the downstream influences with 
change in shape as distinct from that of pure translation. It has been done 
before, and justified empirically or non-dimensionally in several trial-and- 
error or convenient solutions. But in words Kohler seems to be the first who 
has said lag is a function of inflow and attenuation is a function of outflow and 
who, in dimensions, has put one on one side of a line measured one direction 
and the other on the other side of the line measured the other direction — as 
he has in his template of figure 7. 

The writer was not prepared to accept this at face value as rational. But 
being already committed to an empirical usage of that principle, and having ~ 
the means of experimentation at his disposal, he caused two flood waves to be 
generated by purely upstream forces, within the banks of a known stream, and 
without the lateral inflow which confuses natural flood records. Nor was 
there the lack of knowledge of the exact form of the upstream wave as is the 
case with most rain-generated waves. The stream is steep enough so that 
stage-discharge relations are as constant as they can be in ordinary 
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streamflow hydraulics, and the friction factors are so much larger than 
momentum forces and velocity heads that influence of these things on stages 
and discharge is probably much less than it is in the streams for which flood 
routing techniques would be applied. The results are summarized here as 
support for the author’s view (that where the Muskingum routing equations 
apply, i.e. constant time elements and linear storage-discharge relations, the 
effects of flood routing can be suitably approximated by transposition as a 
function of inflow and storage attenuation as a function of outflow). (The 
stream does not have the form or dimensions which would lead to intensifica- 
tion of wave form, so these forms are neither confirmed nor denied by this 
experiment.) 

Two premises are common to the selection of the experimental wave forms, 
the author’s lag routing templates, and the writer’s concept of hydrograph 
analysis and synthesis by basin shaped instantaneous unit-graphs. These two 
premises are that, in a stream where time elements of flow are approximately 
equal regardless of flow magnitude, 


1) Storage effects may be approximated by two kinds of storage in suitable 
combination which tend to attenuate flow and which do not, and, 


2) Duration of flow and increasing and decreasing flows may be approxi- 
mated by small instantaneous impulses of positive or negative change, and 
any hydrograph no matter how complex can be approximated by some se- 
quence of impulses. 


In the writer’s view, premise (1) is thought of first as amounts of storage 
which translates and attenuates, or reductive storage, and some which trans- 
lates and accentuates peaks, with the latter apparently a function of inflow or 
upstream factors and the former a function of outflow or downstream factors. 
This theory is the wedge-storage concept familiar in the Burns-Harkness- 
McCarthy equations. However, these wedges, or two components, can like- 
wise be combined in combinations of two prisms so that tendancies to build 
up and to flatten are balanced and only translation results, and another com- 
ponent providing the net attenuation which is common in most streams (or the 
build up which is possible in some cases and would apparently involve nega- 
tive values of the K-item). This prism combination is Mr. Kohler’s division 
into lag-and-storage values. (The writer would caution that both items repre- 
sent storage, however, and that translatory storage and reductive storage are 
suitable names which retain the concept of physical storage space required in 

-a channel to produce the effects involved). 

With the above premises in mind, the experiments here reported were de- 
signed to release two rectangular wave forms of different magnitude by sud- 
denly opening and suddenly closing gates with suitably sustained periods be- 
tween each change to permit flow to become completely stable and steady 
between changes and steady throughout the zone of the experiment. 
© The resulting wave was observed as river stage at two points, 0.3 and 6.7 
“miles downstream and at five-minute intervals. Using the two stages of posi- 
‘tive steady flow observed and the initial stage as zero flow, an approximate 

rating curve was prepared for each observation point and the stage observa- 
‘tion converted to approximate steady flow in much the manner which is conven- 
tional in stream gaging on the steep streams to which these premises can 


apply. 
| The resulting interpretation is shown on figure 1 hereof. The lag and 
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| reductive storage values are dimensioned only for the downstream point. The 
K values for the descending hydrographs, defined as 


I- 0 
d(I-0) 
dt 


in the general case, or when inflow has ceased as 


(0) 

do 

dt 

are the conventional ones. In this case, unlike natural hydrographs in which 
inflow is never steady except at zero, a determination is also possible on the 
ascending leg, and the values are shown. 

Within the range of many practical uses, the values rising and falling and 
the associated lags are constant. Whether the small differences in the time 
elements as printed on the figure are factual or the faults of methods of ob- 
servation or deduction is not known, and the opportunity to repeat the experi- 
ments suitably has not since presented itself. 

The measure of the whole storage utilized in wave passage is the sum of 
the associated lag and recessional values, successively 258, 298, 310 and 346 
minutes. As volumes these are respectively 0.36, 0.41, 0.43, and 0.48 acre- 
feet per cfs of flow, or nearly 800 acre-feet in the prism of changing water 
level in the range of these observations. The writer wishes to emphasize 
that these are physical volumes, and changes in channel alignment or confine- 
ment could alter them in proportion to those changes. 

This is a stream where the proportion of storage along the channel to flow 
in the channel is about the same throughout the range of in-channel flow rates, 
as has been confirmed in backwater computations using channel cross- 
sections. Hence, this is a stream to which constant time elements apply es- 
Sentially within the channel and the range of these observations. The writer 
cautions against too rapid extension of this experience into streams where 
constant time elements do not apply, and where linear storage-discharge re- 
lations are not a fact. All of the “black marks” on the record of the McCarthy 
flood routing method and on unit hydrograph theory have been put there by 
people who did not first find these constant time elements or linear relation- 
ships to be facts as applicable in their subject streams. Where these rela- 
tionships are factual, the theories work very well. However, where the ap- 
plicable facts are that the storage-discharge curves are wide and complex 
loops, and there is no fixed relationship between stage and flow rate in the 
first place, more basic and more fundamental methods must be used. The 
variable lag-variable reductive storage templates illustrated by Mr. Kohler 
offer more promise than anything else yet available, but, like many other 
loop relationships, may not be constant from flood to flood and wholly precise 
except in the flood from which the loop or template is derived. Nevertheless, 
since the streams to which such loops apply have other severe complications 
of measuring discharge and of shifting relationships, this is only a small 
hazard in the total problem. : 

On the positive side from this cautionary expression, the variable templates 
are capable of indicating the earlier arrival of very large flood peaks and the 
almost complete flattening of small floods just beyond the stages of bank over- 
flow, which phenomena are the characteristics of many larger streams and 
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probably of all aggrading, silt-carrying streams (like the Verdigris River 
which Mr. Kohler used as an example for his template of figure 7 and illus- 
tration of the presented idea). Because the template is adapted to the simplest 
of graphical methods of solution, and involves no reading, tabulation, replotting 
or reprinting, but produces a hydrograph from a hydrograph, it is the most 
promising presently available answer to some difficult flood-time problems. 
One of these is the solution for many streams almost simultaneously, as in 
the emergency of flood-time, and where constantly changing data are being 
received and there is constant demand for the changes in predictions or 
operations expected from the last few hours of rainfall, the breach of some 
levee, or the change in gate setting at some flood control or power dam. 
Under these conditions the efforts of several men must be quickly applied 
under emergency conditions and can be safely done with the templates care- 
fully prepared under the more leisurely conditions of other times. 

The template is the newest and most unique of recent flood tools. It is 
destined to rank for several decades with Sherman’s unithydrograph as a 
basic simplification of essential facts for practical daily work. 
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NORTHEASTERN FLOODS OF 1955: FLOOD CONTROL HYDROLOGY 


Discussion by Leo R. Beard 


LEO R. BEARD.—This paper by Mr. Childs admirably documents the 
hydrologic aspects of a record-breaking flood event and their effects on the 
design of flood-control structures. The writer is particularly interested in 
the author’s discussion relative to flood frequencies. The author justifiably 
feels that “flood frequencies for stations with 100 years of record will give 
fairly reliable frequencies up to 50 years. Beyond this point the reliability 
drops rapidly and beyond 100 years the frequency is fictitious.” However, 
estimates of extreme flood frequencies or probabilities are necessary to an 
adequate economic evaluation of project justification, and proper planning 
therefore requires that satisfactory methods of extrapolating frequency 
curves be developed if at all feasible. The writer has briefly investigated 
this problem along lines explored by Messrs. Russell Morgan, Marshall 
Iakisch and others of the Philadelphia District, Corps of Engineers, in their 
studies of the Delaware River Basin. It is hoped that the following brief 
description of this investigation will shed some light on the evaluation of fu- 
ture frequencies of extreme hurricane floods such as those of August 1955. 

Experience in computing flood frequencies for the New England area has 
been particularly discouraging. While in most of the United States annual 
maximum flood events are apparently distributed in reasonable accord with 
the normal probability curve, there has been strong evidence that the New 
England area experiences an unusually large number of major floods (that is, 
that the frequency curves are skewed in the positive direction). Until the 
1955 floods occurred, there was some debate as to whether the observed ex- 
cess of large floods has been accidental (due to chance) or whether it derives 
from some physical characteristic peculiar to the New England area. With 
the added evidence from the 1955 floods, little doubt remains that some factor 
is operating in New England floods that is not characteristic of floods in the 
rest of the nation. 

The author suggests that “large floods and the smaller events are not re- 
lated and should not be compiled in the same array and analyzed by statistical 
methods,” and “analysis of floods should be divided into two categories: one 
for moderate annual floods, and the second for the larger floods.” The writer 
believes, however, that those factors causing a major event can also combine 
such as to cause a small event, and consequently that some small events 
must be related to the large floods. It is doubtful that there is an abrupt dis- 
continuity in channel conveyance factors such ie they are distinctly differ - 
ent for “large” and “small” floods. 

However, it is agreed that events which are not related (that is, not caused 
by the same general factors) should not be combined in a single frequency ar- 
ray. Since most of the New England floods are caused by extratropical 


a. Proc. Paper 1663, June, 1958, by Elliot F. Childs. 
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cyclones but a few are caused by tropical cyclones (hurricanes), it has been 
suggested that events be segregated into two groups on this basis. In this 
way, perhaps theoretical distribution functions used successfully in other re- 
gions might be applied in the New England area. Such functions are requisite 
to accurate evaluation of extreme flood frequencies. 

This idea which was developed in connection with Delaware River studies, 
was tested by the writer in a study conducted under one of the Civil Works 
Investigations Projects of the Corps of Engineers. Twelve long-record 
stream gaging stations distributed throughout New England were selected for 
the study. The period of record considered was limited to the last 40 years, 
because data on hurricanes during earlier years are not reliable. Maximum 
mean-daily flows were selected in three series as follows: 


a. Annual maximum non-hurricane floods 
b. Annual maximum hurricane floods 
c. Annual maximum floods, regardless of type 


The list of hurricane floods used was compiled from various hurricane re- 
ports. The various reports available were not consistent, and it is felt that, 
particularly in the early years, some of the actual hurricane events have not 
been listed. Hurricane floods used in the study are as follows: 


Hurricane Flood Dates by 10-year Periods 


1916-25 1926-35 1936-45 1946-55 
23 Jul 1916 5 Oct 1927 19 Sep 1936 30 Aug 1949 
1 Oct 1920 11 Aug 1928 21 Sep 1938 20 Aug 1950 
25 Oct 1923 3 Oct 1929 22 Aug 1939 12 Sep 1959 
27 Aug 1924 17 Sep 1932 2 Oct 1939 2 Sep 1952 
1 Oct 1924 24 Aug 1933 2 Sep 1940 15 Aug 1953 
18 Sep 1933 19 Sep 1940 8 Sep 1953 
20 Jun 1934 17 Oct 1943 31 Aug 1954 
9 Sep 1934 2 Aug 1944 12 Sep 1954 
5 Sep 1935 15 Sep 1944 16 Oct 1954 
22 Oct 1944 14 Aug 1955 
27 Jun 1945 18 Aug 1955 
19 Sep 1945 21 Sep 1955 


Since the frequency curve of hurricane floods is based on fewer events 
than the years of record, plotting of the curve was based on the percentage of 
floods rather than on the percentage of years. In order to convert the fre- 
quency in percentage of floods to the frequency per hundred years, the indi- 
cated exceedence frequency must be multiplied by the ratio of the number of 
floods used to the number of years of record. This adjustment of the hurri- 
cane frequency curves is illustrated in Figure 1. 

Frequency curves for one of the stations studied are shown in Figure 1. It 
was found that, in the case of hurricane floods, the plotted data fit the com- 
puted log-normal curves reasonably well in most cases, but occasionally sug- 
gest an upward curvature of the frequency curve. The plotted data of the non. 
hurricane floods almost consistently show an upward curvature of the 
frequency data relative to the log-normal distribution. 
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In order to summarize the data for all twelve stations in such a way as tc 
demonstrate the average curvature, the floods at all stations that were in ex 
cess of the indicated 100-year value were counted. Similarly, the floods tha 
were between the 50-year and 100-year indicated values were counted, and 
forth until the entire range of frequency was covered. Comparison of these 
aggregate numbers with those indicated by the straight line frequency curve 
is shown in Figure 2. It will be noted in Figure 2 that the hurricane-type 
floods evidence a slight upward curvature throughout the range of the fre- 
quency curve. The non-hurricane floods evidence a sharp upward curvature 
only in the upper range of floods. Since it was known that the largest non- 
hurricane general flood that has occurred in New England in the past 300 
years is that of 1936, it was felt that perhaps this one flood may have been 
unduly responsible for the upward curvature shown by the non-hurricane 
floods. In order to examine the effect of this one flood, the 1935 flood event 
were subtracted from the frequency count and the aggregate points replotted 
This comparison, also shown in Figure 2, is somewhat better, but there stil 
is a tendency toward upward curvature. 

When floods are segregated according to cause or type, and where the 
separate frequency curves of annual maximum events are plotted as in this 
study, it is often desirable to obtain a frequency curve of annual floods re- 
gardless of type. In order to do this, the following formula derived in the 
Philadelphia District, Corps of Engineers, is employed: 

P3 _ 12h} “ Bl ot PPS 


100 100 100 10,000 


in which: 


Py = the exceedence frequency per hundred years, annual maximum non 
hurricane floods 


Py = exceedence frequency per hundred years, annual maximum hurrica 
floods 


P3 = exceedence frequency per hundred years, annual maximum floods, 


regardless of type 


Combination of frequencies using this formula is illustrated in Figure 1. 

As might be expected from the experience with the separate hurricane a1 
non-hurricane flood frequency curves described above, combination frequen 
curves computed in the above manner for the 12 stations revealed that the c 
served annual maximum events (regardless of type) at practically all statio 
are more frequent in the range of higher floods than is indicated by the com 
puted combination frequency curves. This may be due to the fact that the 
March 1936 and even the November 1927 non-hurricane floods are much mc 
rare events than would be expected to occur in the 40-year period studied. 
Also, the hurricane flood of August 1955 is considered to be much rarer ths 
can be expected in the 40 years studied. On the other hand, it is possible tk 
the log-normal technique employed, although successful in most regions of. 
United States, is not applicable in New England. Consideration was given tc 
eliminating this disparity by flood segregation other than between hurricane 
and non-hurricane types, but it was concluded that other flood types were n 
distinct or sufficiently independent. ee 

If the results of the study described herein were the only evidence avail: 
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for judging the adequacy of the techniques employed, the excess of experience 
large floods would cast strong doubt on their validity. However, the log- 
normal technique has been successfully employed in most other sections of 
the country, and it is only reasonable to expect a serious excess of large 
floods in some regions and a serious shortage of large floods in other re- 
gions. As a matter of fact, there is a feeling among hydrologists familiar 
with the New England region that the floods of 1936, 1955, and even 1927 are 
indeed very rare events, particularly in the light of historical data. It is be- 
lieved, therefore, that the application of the log-normal technique to the fre- 
quency analysis of hurricane floods and non-hurricane floods separately con- 
stitutes a reasonable method for evaluating the future frequency of extreme 
floods in the New England area. The frequency curves derived by this 
method indicate that the 1955 floods were moderate at 9 of the 12 stations 
studied, but that the flood of 18 August (Hurricane Diane) will be exceeded 
only once in 200 to 300 years at the remaining three stations. 
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